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An  experimental  facility  for  studying  soot  formation  in  high  temperature,  fuel-rich,  laminar,  premixed  flames 
has  been  constructed.  Diagnostics  included  laser  absorption,  thermocouple  particle  densitometry,  and 
thermophoretic  soot  sampling  with  analysis  by  transmission  electron  microscopy.  Single  particles  with 
diameters  as  small  as  3-5  nanometers  were  observed.  Larger  particles  (20-25  nanometers)  agglomerated  to  form 
arge  c  usters.  coflow,  axisymmetric,  laminar  ethylene  diffusion  flame  has  been  studied,  both  experimentally 
mi  computationally.  A  lifted  flame  has  been  selected  to  eliminate  possible  uncertainties  caused  by  the  burner 
lip.  A  two-dimensional,  detailed  soot  growth  model  in  which  the  equations  for  particle  production  are  coupled 
to  e  ow  and  gaseous  species  conservation  equations  has  been  used  to  investigate  soot  production  in  the 
ame.  Detailed  transport  and  finite  rate  chemistry  in  the  gas  phase  was  coupled  with  the  particle  aerosol 
equations  in  the  sectional  representation.  In  comparison  to  measured  data  obtained  using  intrusive  and  non- 
mtrusive  diagnostics,  the  model  predicted  temperature,  flame  height,  and  major  species  very  well.  Peak  benzene 
concentrations  and  soot  volume  fraction  were  predicted  to  within  20%  of  the  experimental  value.  The  predicted 
distribution  of  benzene  was  excellent,  but  the  soot  was  underpredicted  along  the  centerline.  This  deficit  was 
attributed  to  limitations  in  the  PAH  growth  model.  Oxidation  of  particulates  was  dominated  by  reactions  with 
hydroxyl  radical  at  superequilibrium  levels.  Radiation  losses  significantly  effected  predicted  temperatures. 
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FUNDAMENTALS  OF  SOOT  FORMATION 
IN  GAS  TURBINE  COMBUSTORS 


Annual  Progress  Report 
November  15, 1997  to  August  15, 1998 

AFOSR  Contract  F49620-98-C-0008 


I.  Obiectives 

The  overall  objectives  of  this  work  are  to  obtain  necessary  fundamental  data  and  to 
enhance  and  then  validate  modeling  procedures  in  order  to  support  modeling  of  soot 
production  in  practical  gas  turbine  combustors.  Several  focused  tasks  of  this  effort 
include  (1)  extending  measurements  of  surface  growth  rate  constants  to  high 
temperatures,  characteristic  of  rich  zones  in  advanced  combustors;  (2)  modeling  of 
steady,  laminar  diffusion  flames  with  different  fuels  to  assist  in  validating  the  models;  (3) 
advancing  the  modeling  capabilities  for  describing  known  physical  processes  involved  in 
soot  production,  including  carbonization,  ageing,  and  aggregate  formation  to  enable  more 
reliable  extrapolation  of  existing  models;  and  (4)  modeling  of  sooting,  transient  flames  to 
offer  some  physical  understanding  of  processes  controlling  soot  formation  and 
destruction  in  turbulent,  diffusion  flames. 

II.  Status  of  Effort 

The  contract  initiated  at  the  end  of  last  year.  One  major  success  was  a  collaborative  effort 
with  experimentalists  at  Yale,  resolving  a  problem  in  modeling  flame  heights  and  writing 
a  paper  on  modeling  a  coflow  laminar,  ethylene  flame  which  was  accepted  and  presented 
at  the  27th  International  Combustion  Symposium  (see  Appendix  A).  This  work  has  been 
useful  in  identifying  strengths  and  weaknesses  in  both  experiments  and  in  our  model.  A 
new  laboratory  facility  for  studying  soot  formation  in  atmospheric  pressure,  premixed 
flames  has  been  constructed.  Appropriate  diagnostic  tools  have  been  acquired  and  tested 
in  preparation  for  measuring  surface  growth  rates  of  soot  particulates  at  high 
temperatures.  New  procedures  are  being  developed  to  treat  agglomerate  structures  in 
flames.  In  addition,  refinements  to  a  previous  manuscript  have  been  completed  and  it  has 
been  accepted  for  publication  (Appendix  B). 
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TIT.  Accomplishments/New  Findings 

III.  A.  Measurements  of  Surface  Growth  Rates  at  High  Temperatures 

An  important  parameter  for  the  prediction  of  soot  production  is  the  surface  growth  rate. 
Virtually  all  reliable  measurements  of  this  parameter  have  been  obtained  over  the 
temperature  range  of  1500-1950K.  Many  authors  simulate  this  rate  using  a  rate  constant 
which  increases  exponentially  with  temperature  (Arrhenius  formulation).  This 
assumption  is  acceptable  under  standard  conditions  since  soot  formation  usually  occurs  in 
diffusion  flames  at  temperatures  below  1800K.  Ref  1  and  2  have  separately  developed 
temperature-dependent  expressions  showing  a  substantial  fall-off  in  the  net  rate  constant 
for  surface  growth  above  about  1850K.  These  formulations  were  in  part  developed  based 
on  limited  data  at  elevated  temperatures  (see  Ref  3).  Reliable  data  at  temperatures  above 
2000K  is  very  limited  and  hence  extrapolation  of  this  fall-off  trend  is  highly  uncertain. 
This  problem  is  a  concern  since  for  advanced,  high  performance  gas  turbine  combustors 
the  ‘sooting’  rich  zone  will  be  at  temperatures  well  in  excess  of  2000K.  This  problem  is 
depicted  in  Figure  1  in  which  equilibrium  temperatures  for  several  different  conditions 
are  plotted.  The  temperatures  computed  with  the  heated  ethylene  fuel  and  high  air  preheat 
can  be  taken  as  advanced  combustor  conditions.  Since  most  data  is  obtained  for  the  lower 
temperature  zone,  there  is  an  obvious  concern  in  extrapolating  existing  models  to  the 
conditions  of  the  advanced  combustors. 


Figure  1 


Comparison  of  Equilibrium  Flame 
Temperatures 
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To  provide  the  needed  data,  a  series  of  fuel-rich,  laminar  premixed  flames  are  being 
examined,  analogous  to  several  previous  studies  (see  for  example,  Ref.  4-5).  Soot 
evolution  in  ethylene/oxygen/argon  flames  with  C/O  ratios  from  0.56  to  0.92  have  been 
characterized  using  laser  extinction,  thermocouple  particle  densitometry  (TPD),  and 
thermophoretic  sampling  followed  by  transmission  electron  microscopy  (TEM). 

Initial  measurements  of  soot  volume  fraction  using  TPD  were  in  good  agreement  with 
those  made  using  laser  extinction.  The  measurements  were  made  on  a  custom  hastalloy 
burner  containing  513  x  1.2mm  diameter  holes  spaced  2.0mm  o.c.  (28%  open  area)  which 
is  surrounded  by  a  co-flow.  Concern  over  the  uniformity  and  ultimate  flame  holding 
ability  at  future  flow  conditions  of  this  burner  design  prompted  a  redesign  to  a  sintered 
metal  surface  burner.  Results  from  the  new  burner  will  be  used  to  finalize  our 
implementation  of  the  TPD  technique  which  depends  on  proper  analysis  of  the  soot 
deposition  stages  of  the  thermocouple  junction. 

Rapid  insertion  and  withdrawal  of  3mm  diameter  carbon-coated  TEM  grids  accomplished 
thermophoretic  sampling  in  the  flame.  Subsequent  inspection  of  the  deposits  using  TEM 
indicated  the  presence  of  individual  (as  small  as  5  nm  diameter)  and  agglomerated 
particles.  The  agglomerate  size  increased  with  downstream  distance  while  (consistent 
with  prior  work)  the  primary  particle  size  in  the  agglomerates  remained  about  20nm.  A 
computer  controlled  sampling  system  is  nearing  completion  and  further  refinement  of  the 
TEM  analysis  is  required. 


III.  B.  Modeling  of  Soot  Formation  in  Diffusion  Flames 

In  addition  to  work  in  modeling  opposed  jet  diffusion  flames  (see  Ref.  6-7),  significant 
accomplishments  in  the  modeling  of  a  coflow,  laminar  diffusion  flame  model  have  been 
achieved.  Most  studies  in  which  detailed  chemical  kinetics  are  coupled  with  detailed  soot 
models  have  been  focused  on  one-dimensional  problems.  Ref.  8  and  9  have  modeled 
laminar  jet  diffusion  flames  using  monodisperse  soot  formation  models  with  skeletal 
reaction  mechanisms.  In  this  program,  we  have  incorporated  the  sectional  soot  formation 
and  radiation  models  developed  in  Ref.  1  and  6  into  the  laminar,  axisymmetric,  diffusion 
flame  code  (see  Ref.  10)  for  a  cylindrical  fuel  stream  surrounded  by  a  coflowing  oxidizer 
jet.  This  recent  work  is  described  in  detail  in  Ref.  11  (Appendix  B)  for  an  attached 
methane  flame  and  in  Ref.  12  (Appendix  A)  for  a  diluted,  lifted  ethylene  flame.  For 
methane,  a  modified  GRIMech  2.11,  with  all  NO  related  reactions  deleted  and  some 
benzene  formation  and  destructions  steps  added,  was  used  for  the  computations.  The 
ethylene  mechanism  was  derived  from  GRIMech,  based  on  comparisons  to  PSR  data  and 
ignition  delay  times.  We  have  collaborated  with  colleagues  at  Yale  who  have  been 
investigating  these  flames  using  a  variety  of  experimental  techniques  (thermocouples, 
TPD,  quartz  probes  with  on-line  mass  spectrometry  and,  in  the  case  of  the  ethylene  study, 
planar  laser  imaging  using  Rayleigh  scattering  and  laser-induced  incandescence). 
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For  the  attached  methane  flame,  qualitative  agreements  were  good,  but  several 
quantitative  differences  between  the  model  and  experiment  were  apparent.  Most  of  these 
were  attributable  to  an  inability  to  model  accurately  the  bulk  flow  features  of  the  flame 
(particularly  flame  height  was  over  predicted  and  peak  flame  temperatures  were  under 
predicted.  Given  these  important  differences,  the  agreement  of  species  concentrations  and 
soot  profiles  were  quite  good.  However,  quantitative  comparisons  between  the  measured 
and  predicted  soot  profiles  could  not  be  made  (predictions  were  about  a  factor  of  four 
low).  Predicted  distributions  of  particle  size  as  a  function  of  radius  is  shown  in  Fig.  2  at  a 
height  of  3  cm  (approximately  one  cm  below  the  height  at  which  the  peak  centerline 
value  is  obtained).  These  distributions  can  be  shown  to  be  representative  of  the  two  soot 
producing  regions  of  the  flame.  One  along  the  wings  (just  inside  of  the  flame  front)  where 
peak  volume  fractions  are  attained  and  the  other  along  the  centerline  where  much  slower 
growth  and  particle  sizes  are  observed. 


Figure  2.  Contributions  of  various  soot  mass  classes  to  the  total  soot 
volume  fraction  as  a  function  of  radius  at  a  height  of  3  cm.  Labels  for  the 
mass  classes  include  the  minimum  and  maximum  diameters  (nm). 

The  diluted  lifted  ethylene  flame  was  selected  in  order  to  avoid  possible  complications  of 
interactions  between  the  burner  lip  and  the  flame,  a  possible  source  of  the  above 
discrepancies.  In  addition,  optical  diagnostics  were  included  to  add  additional 
information  to  help  resolve  any  differences  between  the  model  and  physical  probing  of 
the  flame.  Overall,  the  agreements  between  the  model  and  the  experiments  were 
extremely  good  for  many  features  of  the  flame.  Between  the  three  solutions  for  the 
ethylene  flame  (one  model  and  two  experiments),  usually  at  least  two  of  the  methods 
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agreed  well  (both  spatial  distributions  and  quantitative  values).  A  comparison  of  the  gas 
sampling  and  model  predictions  for  benzene  is  shown  in  Fig.  4  of  Appendix  A,  with 
identical  scales  for  both  contour  plots.  Very  good  qualitative  and  quantitative  agreement 
is  obtained.  A  limitation  in  the  soot  model  can  be  illustrated  by  examining  Fig.  3a,  b,  and 
c,  which  compare  contour  plots  of  total  soot  volume  fraction  for  the  model  and  the  two 
experimental  methods.  While  the  maximum  values  of  soot  volume  are  reproduced  well, 
the  predicted  soot  profiles  peak  in  the  wing  tips,  while  both  experimental  methods  (TPD 
and  LII)  peak  along  the  centerline.  A  preliminary  analysis  of  this  failure  in  the  model  is 
its  inability  to  reproduce  accurately  the  processes  of  early  PAH  formation  occurring  along 
the  centerline  of  the  burner.  Additional  work  in  this  area  is  proceeding  as  are  initial 
efforts  to  model  the  ‘Santoro’  ethylene  diffusion  flame. 


Figure  3.  Comparison  of  predicted  and  measured  soot  volume  fractions 


r  (cm)  r(cm)  r(cm) 


III.  C .  Modeling  the  Growth  and  Oxidation  of  Soot  Aggregates 

In  modeling  the  growth  of  soot  particles,  it  has  been  conventional  to  assume,  for 
simplicity,  that  the  particles  consist  of  isolated  spheroids,  which  undergo  surface  growth 
and  oxidation,  and  coalesce  with  one  another  on  impact.  The  experimental  evidence, 
however,  is  that  while  the  coalescing  spheroid  picture  is  accurate  for  very  small  (young) 
particles,  older  soot  particles  tend  to  consist  of  chainlike  aggregates  of  primary  spheroids 
which  fuse  together  on  impact  but  do  not  coalesce.  The  sectional  aerosol  dynamics 
representation  that  has  been  used  to  model  the  growth  of  spheroids  is  being  adapted  to 
model  aggregates,  as  well. 
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The  algorithm  being  adapted  for  use  in  soot  kinetics  has  a  provision  for  a  discrete  size 
range  of  precursor  particles  whose  masses  are  multiples  of  some  monomer  growth 
species.  The  discrete  size  range  is  joined  to  an  isolated  spheroid  size  range  whose 
boundaries  vary  linearly  on  a  logarithmic  scale;  these  “liquid”  particles  are  assumed  to 
coalesce  on  collision.  Beyond  an  arbitrarily  determined  size,  particles  are  assumed  to 
fuse  with  one  another,  leading  to  aggregate  formation.  The  conventional  spheroid 
sectional  analysis  tracks  the  total  mass  density  within  set  sectional  boundaries;  the 
extended  sectional  analysis  tracks  not  only  the  mass  density,  but  also  the  number  of 
primary  particles  within  a  section.  Aggregate-aggregate  collisons  leading  to  larger 
aggregates  are  prescribed  at  rates  based  on  an  assumed  fractal  dimension  for  the  clusters. 
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V.  Publications 


C.  S.  McEnally,  A.  M.  Schaffer,  M.  B.  Long,  L.  D.  Pfefferle,  M.  D.  Smooke,  M.  B. 
Colket,  and  R.  J.  Hall,  “Computational  and  Experimental  Study  of  Soot  Formation  in  a 
Coflow,  Laminar  Ethylene  Diffusion  Flame,”  To  be  published  in  the  Twenty-Seventh 
Symposium  (International)  on  Combustion  (1999)  (see  Appendix  A.) 

M.  D.  Smooke,  C.  S.  McEnally,  L.  D.  Pfefferle,  R.  J.  Hall,  and  M.  B.  Colket, 
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VI.  Interactions/Transitions 
VI.  A.  Meetings 

On  October  27-29,  1997,  M.  Colket,  D.  Liscinsky,  R.  Hall  and  M.  D.  Smooke  attended 
the  Fall  Technical  Meeting  of  the  Eastern  States  Section  of  the  Combustion  Institute  at 
Hartford,  CT.  A  paper  entitled  “Interpretations  of  a  Computational  Study  of  Soot 
Formation  Coflow,  Laminar  Diffusion  Flame”  by  M.  Smooke,  R.  Hall,  and  M.  Colket 
was  presented. 

On  March  29- April  1,  1998,  M.  Colket  attended  the  20th  Annual  Combustion  Research 
Conference  of  the  US.  DOE,  Basic  Energy  Sciences,  Dallas,  TX.  This  meeting  was  held 
in  conjunction  with  the  American  Chemical  Society  National  Meeting. 

On  April  21st’  1998,  M.  Colket  went  to  NASA  LeRC  and  led  discussions  on  particulate 
measurements  and  soot  modeling. 

On  August  2-7, 1998,  M.  Colket,  R.  Hall,  and  M.  Smooke  attended  the  27th  International 
Symposium  on  Combustion  in  Boulder,  CO.  A  paper  entitled  “Computational  and 
Experimental  Study  of  Soot  Formation  in  a  Coflow,  Laminar  Ethylene  Diffusion  Flame” 
was  presented.  In  addition,  a  (related,  but  separately  supported)  poster  paper  entitled 
“Computational  and  Experimental  Study  of  Cj  to  C6  Hydrocarbons  in  a  Steady, 
Axisymmetric  Flame”  was  presented. 

VI.  B.  Advisory  Functions 

None 

VI.  C.  Transitions 

In  1997  and  1998,  M.  Colket  refined  software  for  rapidly  calculating  equilibrium 
temperatures  and  NOx  formation  and  transferred  it  to  Dr.  Saadat  Syed  ((561)796-3560)  of 
Pratt  and  Whitney  for  imbedding  into  a  CFD  design  code. 
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VII.  Record  of  Inventions 

There  were  no  inventions  during  this  reporting  period.  (DD  Form  882,  Interim  Patent 
Report  to  be  submitted  under  separate  cover) 

VIII.  Honors/Awards 

M.  Colket,  as  a  member  of  a  CFD  development  team,  has  been  nominated  and  is  a 
finalist  for  the  1998  Pratt  and  Whitney  Leadership  Award. 
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Coflow,  Laminar  Ethylene  Diffusion  Flame 


A-l 


COMPUTATIONAL  AND  EXPERIMENTAL  STUDY  OF 
SOOT  FORMATION  IN  A  COFLOW,  LAMINAR  DIFFUSION  FLAME 


M.  D.  Smooke,  C.  S.  McEnally,  L.  D.  Pfefferle 
Yale  Center  for  Combustion  Studies 
Yale  University,  New  Haven,  CT  06520 
and 

R.  J.  Hall  and  M.  B.  Colket 
United  Technologies  Research  Center 
E.  Hartford,  CT  06108 


Abstract 

A  detailed  soot  growth  model  in  which  the  equations  for  particle  production  have  been 
coupled  to  the  flow  and  gaseous  species  conservation  equations  has  been  developed  for  an 
axisymmetric,  laminar,  coflow  diffusion  flame.  Results  from  the  model  have  been  com¬ 
pared  to  experimental  data  for  a  confined  methane-air  flame.  The  two-dimensional  system 
couples  detailed  transport  and  finite  rate  chemistry  in  the  gas  phase  with  the  aerosol  equa¬ 
tions  in  the  sectional  representation.  The  formulation  includes  detailed  treatment  of  the 
transport,  inception,  surface  growth,  oxidation,  and  coalescence  of  soot  particulates.  Ef¬ 
fects  of  thermal  radiation  and  particle  scrubbing  of  gas  phase  growth  and  oxidation  species 
are  also  included.  Predictions  and  measurements  of  temperature,  soot  volume  fractions 
and  selected  species  are  compared  over  a  range  of  heights  and  as  a  function  of  radius. 
Flame  heights  are  somewhat  overpredicted  and  local  temperatures  and  volume  fractions 
are  underpredicted.  We  believe  the  inability  to  reproduce  accurately  bulk  flame  parameters 
directly  inhibits  the  ability  to  predict  soot  volume  fractions  and  these  differences  are  likely 
a  result  of  uncertainties  in  the  experimental  inlet  conditions.  Predictions  of  the  distribu¬ 
tions  of  particle  sizes  indicate  the  existence  of  (relatively)  low  molecular  weight  species 
along  the  centerline  of  the  burner  and  trace  amounts  of  the  particles  which  escape  from 
the  flame,  unoxidized.  Oxidation  of  particulates  is  dominated  by  reactions  with  hydroxyl 
radicals  which  attain  levels  approximately  ten  times  higher  than  calculated  equilibrium 
levels.  Gas  cooling  effects  due  to  radiative  loss  are  shown  to  have  a  very  significant  effect 
on  predicted  soot  concentrations. 
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COMPUTATIONAL  AND  EXPERIMENTAL  STUDY  OF 
SOOT  FORMATION  IN  A  COFLOW,  LAMINAR  DIFFUSION  FLAME 


1.  Introduction 

In  the  last  few  years,  we  have  witnessed  a  shift  in  combustion  research.  The  push  for 
higher  combustion  efficiency  in  propulsion  applications  that  dominated  much  of  research 
in  the  past  few  decades  is  gradually  being  replaced  by  a  drive  towards  cleaner  combustion. 
This  is  a  direct  result  of  environmental  consciousness  and  it  has  been  translated  into  stricter 
air  quality  legislation.  Although  it  originated  as  a  reaction  to  regulatory  pressure,  research 
in  the  area  of  pollutant  formation  and  control  will  become  economically  indispensable  to 
the  export  of  combustion-related  technologies  and  products  worldwide.  In  particular,  as 
emissions  legislation  becomes  more  restrictive,  a  detailed  understanding  of  soot  formation 
and  unburned  hydrocarbons  in  flames  will  be  critical  for  the  design  of  pollutant  abatement 
strategies  and  for  the  preservation  of  the  competitiveness  of  combustion  related  industries. 

Modeling  soot  formation  in  practical  combustion  systems  is  an  extremely  challenging 
problem.  Reacting  flow  problems  for  simple  fuels  such  as  methane  can  require  the  solution 
of  more  than  50  chemical  species  in  addition  to  the  temperature  and  the  fluid  dynamic 
variables.  The  inclusion  of  soot  inception,  growth  and  oxidation  processes  in  these  models 
greatly  increases  the  level  of  complexity  of  the  problem.  In  fact,  due  to  the  number  of 
dependent  unknowns  and  the  resolution  with  which  reaction  fronts  need  to  be  resolved, 
gas  phase  combustion  calculations  with  detailed  chemistry  have  only  recently  moved  from 
one-  to  two-  and  some  three-dimensional  systems.  The  increased  effort  with  which  mul¬ 
tidimensional  combustion  phenomena  have  been  studied  numerically  is  attributable  in 
part  to  advances  in  the  development  of  computational  algorithms,  and  in  part  to  higher 
speed/larger  memory  workstations.  This  combination  has  enabled  the  combustion  scien¬ 
tist  to  investigate  chemically  reacting  systems  that  were  computationally  infeasible  only  a 
few  years  ago.  Nevertheless,  the  investigation  of  soot  formation  with  detailed  chemistry 
in  a  generic  multidimensional  configuration  is  still  beyond  our  current  computational  abil¬ 
ity.  The  laminar  diffusion  flame,  however,  provides  a  natural  environment  in  which  one 
can  investigate  the  interaction  of  soot  formation  with  detailed  gas-phase  chemistry  in  a 
multidimensional  system. 

Most  studies  using  detailed  chemical  kinetics  and  coupled  models  of  soot  production 
and  oxidation  have  focused  on  one-dimensional  geometries.  A  few  recent  studies  [Kaplan 
et  a L,  1996,  Kennedy  et  a/.,  1996]  have  focused  on  jet  diffusion  flames,  using  simplified, 
monodisperse,  soot  formation  models  with  skeletal  kinetic  mechanisms.  In  this  paper  we 
modify  the  sectional  soot  formation  model  developed  in  Colket  and  Hall  (1994)  and  Hall  et 
a 1.  (1997)  for  incorporation  into  a  laminar,  axisymmetric,  diffusion  flame  in  which  a  cylin- 
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drical  fuel  stream  is  surrounded  by  a  coflowing  oxidizer  jet.  Computationally, 
we  employ  a  velocity-vorticity  model  (Ern  et  a L,  1995)  in  which  the  governing  conser¬ 
vation  equations  are  solved  with  detailed  transport  and  finite  rate  chemistry  submodels 
to  predict  the  temperature,  species  mass  fractions  and  velocity  fields  as  functions  of  the 
two  independent  coordinates.  Appropriate  sectional  equations  are  included  for  describing 
the  different  soot  particle  size  classes  with  a  convective,  diffusive  and  chemical  production 
balance.  These  equations,  along  with  surface  growth,  oxidation  and  radiation  (from  both 
gas  and  particulate  species)  are  fully  integrated  with  the  governing  equations.  A  discrete 
solution  is  obtained  on  a  two-dimensional  grid  by  employing  Newton’s  metho-  with  adap¬ 
tive  mesh  refinement.  In  many  models  diffusion  in  the  axial  direction  is  neglected:  here, 
we  consider  the  fully  elliptic  problem.  Experimentally,  we  utilize  thermocouple  particle 
densitometry  and  mass  spectrometry  to  obtain  profiles  of  the  soot  volume  fraction  and  the 
species  concentrations  at  various  heights  in  the  flame  (McEnally  and  Pfefferle,  1997). 

2.  Problem  Formulation 

While  the  ultimate  goal  in  combustion  modeling  is  the  solution  of  three-dimensional 
turbulent  reacting  flows  with  finite  rate  chemistry,  there  are  still  important  less  complex 
systems  that  can  be  analyzed  in  detail  with  current  computational  resources.  The  axisym- 
metric,  laminar,  diffusion  flame  is  one  such  configuration  (Smooke  et  a J.,  1990,  1992,  1996; 
Xu  et  a J.,  1993;  Ern  et  a/.,  1995).  Diffusion  flames  are  typical  of  most  practical  combustion 
devices.  The  ability  to  predict  the  coupled  effects  of  complex  transport  phenomena  with 
detailed  chemical  kinetics  even  in  laminar  systems  is  useful  in  the  modeling  of  turbulent 
reacting  flows  and  in  understanding  the  processes  by  which  pollutants  are  formed. 

2.1  Soot  Modeling 

Soot  kinetics  are  modeled  as  coalescing,  solid  carbon  spheroids  undergoing  surface 
growth  in  the  free  molecule  limit.  The  particle  mass  range  of  interest  is  divided  into  sec¬ 
tions  (Gelbard  and  Seinfeld,  1980),  and  an  equation  is  written  for  each  section  including 
coalescence,  surface  growth,  and  oxidation.  In  the  sectional  representation,  the  sectional 
mass  boundaries  vary  linearly  on  a  logarithmic  scale.  Sectional  analysis  makes  it  possible 
to  obtain  the  particle  size  distribution  without  apriori  assumptions  about  the  form  of  the 
distribution.  For  the  first  bin,  an  inception  source  term  is  included.  The  incorporation 
of  these  equations  into  a  transport/conservation  equation  for  each  section  includes  both 
thermophoresis  and  an  effective  bin  diffusion  rate.  In  the  gas  species  conservation  equa¬ 
tions,  provision  is  made  for  complex  chemistry;  additional  source  terms  are  included  to 
account  for  scrubbing  or  generation  of  gaseous  species  arising  from  the  particle  growth 
and  oxidation  processes.  The  gas  and  soot  are  additionally  coupled  through  non-adiabatic 
radiative  loss  from  both  the  gas  and  the  soot  in  the  optically-thin  approximation.  For  one 
atmosphere  calculations,  the  sectional  coefficients  for  transport,  coalescence  and  surface 
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processes  are  calculated  in  the  free  molecule  limit  (Gelbard  and  Seinfeld,  1980),  with  sur¬ 
face  growth  and  oxidation  proportional  to  particle  surface  area.  In  the  free  molecule  limit, 
the  dependencies  of  the  sectional  coefficients  of  the  particle  mass  in  its  surrounding  gas 
medium  are  factored  in  such  a  way  that  they  only  have  to  be  evaluated  once.  Oxidation 
of  soot  by  O2  and  OH  is  treated  as  described  in  (Hall  et  aJ.,  1997).  The  inception  model 
employed  here  is  based  on  an  estimate  of  the  formation  rate  of  two-  and  three-ringed  aro¬ 
matic  species,  and  is  a  function  of  local  acetylene,  benzene,  phenyl  and  molecular  hydrogen 
concentrations.  Assuming  steady-state  values  of  intermediates,  and  also  that  H2  >>  C2H2, 
the  rate  of  production  of  the  polyaromatic  species  can  be  estimated  to  be  (Hall  et  a I.,  1997) 

d[Cl0H7]  _  10n.Me(— 4378/T) [^l![C6g5]  cc/mole/sec,  (2.1) 

dt  [H2] 

and 


d[ClAHl0]  _  I0i2.5oe(-6390/T)1^2^2]^6^6][c6>f5]  Cc/mole/sec,  (2.2) 

dt  [H2] 

where  the  gas  phase  concentrations  and  temperatures  are  evaluated  at  local  conditions. 
With  the  further  assumptions  that  inception  is  limited  by  the  formation  of  polyaromatics 
and  oxidation/decomposition  of  such  species  can  be  neglected,  the  inception  rate,  5,,  in 
grams/cc/sec,  was  initially  assigned  to 

S,  =  127  x  <f[C'°gj  +  178  x  (2.3) 

dt  dt 

where  the  constants  (molecular  weights)  are  provided  to  convert  from  molar  to  mass  units. 
The  second  term  in  equation  2.3  plays  a  negligible  role  in  methane-air  diffusion  flames  but 
has  been  included  here  for  generality.  The  contributions  from  both  inception  processes  are 
incorporated  in  the  first  sectional  bin,  whose  lower  mass  boundary  is  set  equal  to  the  mass 
of  the  smallest  inception  species  (See  Section  5). 

This  inception  model  contrasts  with  those  utilized  by  Frenklach  and  coworkers,  see 
Markatou,  et  al.,  (1993)  and  references  contained  therein.  Originally,  such  models  included 
detailed  chemical  kinetics  to  describe  the  formation  of  very  large  molecular  weight  species. 
Their  inception  model  has  subsequently  been  modified  to  include  the  dimerization  of  much 
smaller  PAH  species,  such  as  pyrene  (Ci6-ff\o)-  Our  use  of  smaller  species  for  inception  likely 
overestimates  the  initial  rates  of  inception,  but  then  also  overestimates  rates  of  coalescence 
which  reduces  surface  area  and  surface  growth.  While  these  effects  are  not  counterbalanced, 
the  uncertainties  in  PAH  formation  models  and  their  role  in  soot  nucleation  do  not  yet 
justifjr  the  use  of  the  more  detailed  approach. 
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The  surface  growth  rate  is  nominally  that  of  Harris  and  Weiner  (1983),  who  report 
a  rate  proportional  to  acetylene  concentration  near  1650  K.  We  have  imposed  a  31.8 
kcal/mole  activation  energy  as  suggested  by  Hura  and  Glassman  (1988),  and  have  empiri¬ 
cally  multiplied  the  surface  growth  Arrhenius  A-factor  by  a  factor  of  two.  'The  enhancement 
by  a  factor  of  two  brings  the  Harris  and  Weiner  rate,  determined  from  studies  on  premixed 
ethylene  flames,  into  better  agreement  with  results  obtained  in  diffusion  flames  (Axelbaum 
et  ai.,  1988,  Kennedy  et  a i.,  1990,  Sunderland  and  Faeth,  1996),  where  somewhat  larger 
surface  growth  rates  have  been  measured. 

It  is  not  suggested  that  the  soot  growth  model  used  here  represents  a  universal  solu¬ 
tion  to  the  problem  of  soot  formation  because  comparisons  with  soot  growth  from  other 
fuels  have  not  been  performed.  The  possibility  of  aging  of  the  surface  growth  rate  is  not 
accounted  for  here,  nor  is  aggregate  formation.  These  effects  should  be  included  m  future 
modeling  efforts.  Fundamental  questions  remain  regarding  inception,  temperature  depen¬ 
dence  of  surface  growth  and  precise  dependence  of  surface  growth  on  gas  phase  species 
concentrations.  Whether  soot  growth  rates  derived  from  premixed  and  opposed  jet  studies 
are  self-consistent  is  still  an  open  issue  as  well. 

The  number  of  soot  sections  required  for  convergence  of  soot  size/density  parameters 
must  be  examined  in  each  problem,  although  evidence  points  to  the  relative  magnitudes 
of  surface  growth  and  inception  as  determining  this  number.  When  the  contribution  of 
inception  to  total  volume  fraction  is  significant,  although  not  necessarily  dominating,  the 
number  of  sections  required  can  be  very  small  (Hall  et  ai.,  1997),  with  larger  numbers 
necessary  when  the  contribution  of  surface  growth  to  volume  fraction  becomes  very  much 
larger  than  that  from  inception. 

2.2  Governing  Equations 

If  we  combine  the  gas-phase  diffusion  flame  model  which  employs  a  velocity- vorticity 
formulation  (Ern  et  ai.,  1995)  with  the  sectional  approach  presented  in  Hall  et  al.  (1997), 
the  coupled  set  of  governing  partial  differential  equations  can  be  written  in  the  form 

Radial  Velocity 

d2vr  d2vr  _  du  d_  (^dA  _  d_  fvA  _  d_  /V  V\  /2.4) 

~drJ+  dz2  ~  dz  dr\pdz)  dr  \r  J  dr  \  p  dr )  ’ 
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Axial  Velocity 


d'-v,  d'-v,  _  _  _a  /»,  M  _  i.  I'M  _  i.  / JV ^  _ 

ftT+  flri  _  5z  3z  l  p  dz )  dz  \  r  )  dz\pdrj' 
Vorticity  Transport 


d 

dz 


d  fujdip\ 
dz\r  dr  j 


a 

(  u>  dip  \ 

d 

’  3  A  r-Y 

dr 

{-rTz). 

dr 

dr  \  r  /  _ 

-r9^-rv(^)-Mp) 


(2.5) 


(2.6) 


Species 


pVr 


dV \ 

dr 


+  +  ;l :MVt,)  +  £:{rYtv..)- 


Wk(wk+w’k)  =  0, 


dz 

k  =  1,2,...,  A', 


(2.7) 


Sections 


dYk 

pv'~d7 


+  pv^  +  ~{rpYk(yZ  +  Vrr))  +  Q~z  (pYkW,  +  VTt)) 
-qk  =  0,  k  =  K  +  l,K  +  2,...,K  +  M, 


(2.8) 


Energy 


dT 


r-fr +  pv-Tz , 


dT\  i  a 


_  a.  (xst\ 

rdr\Xdr)  dz  \X  dz  )  + 


E Pc>-Y‘  (  j;Vl1  +  T:Vt‘)  +  h‘lVl‘(-’i’t  +">;)-  V  •  9'  =  °’ 


(2.9) 


Equation  of  State 


o=£- (2-10) 

p  RT 
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The  system  is  closed  with  appropriate  boundary  conditions  on  each  side  of  the  computa¬ 
tional  domain.  For  a  confined  coflowing  flame  we  have 

Axis  of  Symmetry  (r  =  0): 


_dv1_du  =  9YjL  =  ?Z  =  0  k  =  1,2, . . . ,  K  +  M, 

Vr  ~  dr  ~  dr  dr  dr 

Exit  (z  — »  oo): 

an=5T=  k  =  l,2,...,K  +  M, 

Vr~  dz  ~  dz  dz  dz 

Inlet  (z  =  0): 
r  <  Ri 

VT  =  0, 

Vz  =  Vf, 

dvr  dvz 

U  =  lh~lh' 

Yk  =  YkF,  + 

T  =  Tf, 

Rj  <  r  <  R0 

Vr  =  0, 

Vt  =  Vox , 

dvT  dvz 

IJ  =  Th~!h’ 

n  =  n„„  k  =  i,2,...,K+M, 

T  =  Tax , 


(2.11) 


(2.12) 


(2.13) 


(2.14) 


Wall  (r  =  Jio): 


Vr  =Vz=  0, 


^  =  0,  k  =  l,2,...,K  +  M, 

dr 


T  =  T. 


wall  • 


(2.15) 


The  subscripts  I  and  0  refer  to  the  inner  jet  and  the  outer  jet,  respectively,  and  vF,vox, 
YkF,Ykox,TF,Tox  and  Twalt  are  specified  quantities. 
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In  addition  to  the  variables  already  defined,  A  represents  the  total  number  of  gas 
phase  species  and  M,  the  number  of  soot  (size)  sections;  T,  denotes  the  temperature;  Yk, 
the  mass  fraction  of  the  *“  gas-phase  species  (k  <  A)  or  the  (k  -  A)  soot  size  class 
Ik  >  A);  p,  the  pressure;  ty  and  v„  the  velocities  of  the  fluid  mixture  jn  the  radial  and 
axial  directions,  respectively;  p,  the  mass  density;  Wk1  the  molecular  weight  of  the  k 
species;  W,  the  mean  molecular  weight  of  the  mixture;  R,  the  universal  gas  constant;  A, 
the  thermal  conductivity  of  the  mixture;  cr,  the  constant  pressure  heat  capacity  of  the 
mixture;  c,„  the  constant  pressure  heat  capacity  of  the  k‘h  species;  u>»,  the  molar  rate 
of  production  of  the  k,h  species  per  unit  volume;  the  molar  rate  of  production  o  e 
k’h  species  per  unit  volume  due  to  scrubbing/replenishment  by  the  soot  growth/oxidation 
processes;  q„  the  rate  of  change  of  section  k  due  to  inception,  surface  growth  oxidation 
and  coalescence;  ht,  the  specific  enthalpy  of  the  V  species;  g,  the  gravitational  constant; 
„  the  viscosity  of  the  mixture;  V„,  and  Vt„  the  diffusion  velocities  of  the  k  species  in 
the  radial  and  axial  directions,  respectively;  VTr,  VT„  V;„  V;„  sectional  thermophoretic  and 
diffusion  velocities  of  the  k'h  soot  size  class  and  V  ■  «,  the  divergence  of  the  net  radiative 
flux  for  gas  bands  and  soot  in  the  optically-thin  limit. 


We  write  the  diffusion  velocities  in  the  r  and  z  directions  in  the  form 

H,  = 

K.  =  -(l/A,)Hi^ 


*  =  1,2,. 

..,/G 

(2.16) 

*  =  1,2,. 

(2.17) 

and  Dk  is 

related  to  the  binary  diffusion 

VV  llUi  — - - -  - 

coefficients  through  the  relation  (see,  e.g.,  Curtiss  and  Hirschfelder,  1949) 


,  (1  ~  y>) 


(2.18) 


The  sectional  thermophoretic  velocities  in  the  free  molecule  regime  are  given  by  (Hall  et 
at.,  1997) 


u  1  dT 

vTr  = 

u  1  dT 

Vtz  =  _0-55pT 


(2.19) 

(2.20) 


The  sectional  diffusion  velocities  are  written  as  in  (2.16-2.17)  with  a  mass  weighted  mean 
diffusion  coefficient  for  each  bin  (Hall  et  a 1.,  1997).  The  binary  diffusion  coefficients,  the 
viscosity,  the  thermal  conductivity  of  the  mixture,  the  chemical  production  rates  and  the 
thermodynamic  quantities  were  evaluated  using  vectorized  and  highly  optimized  transport 
and  chemistry  libraries  (Giovangigli  and  Darabiha,  1987). 

Anticipating  that  radiative  losses  could  have  a  significant  influence  on  soot  levels  com¬ 
pared  to  a  nonradiative  system,  we  have  included  an  optically-thm  radiation  model  m  our 
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calculations  (Hall,  1993,  1994).  We  assume  that  for  methane-air  mixtures  the  only  sigm  - 
cant  radiating  species  are  H,0,  CO,  CO 2  and  soot.  By  utilizing  an  optically-thm  limit  m 
which  self-absorption  of  radiation  is  neglected,  the  divergence  of  the  net  radiative  flux  can 

be  written 

V  -qR  =  CfvT5  +  (2>21) 


where  /,  is  the  soot  volume  fraction  and  Ihtk  is  the  Planck  function  evaluated  at  the 
gas  band  centers  of  the  contributing  vibration-rotation  or  pure  rotational  bands  whose 
integrated  intensities  are  given  by  The  value  of  C  is  taken  to  be  4.243E-10  (giving,  a 
power  density  in  watts/cc  for  T  in  Kelvin),  and  was  derived  from  Grosshandler  (1993). 

3.  Experimental  Methods 

3.1  Burner  and  Reactants 

Atmospheric  pressure,  overventilated,  axisymmetric,  coflowing,  nonpremixed  laminar 
flames  were  generated  with  a  burner  in  which  the  fuel  flows  from  an  uncooled  12  mm  diam¬ 
eter  vertical  brass  tube  and  the  oxidizer  flows  from  the  annular  region  between  this  tu  e 
and  a  102  mm  diameter  concentric  tube  (see  Figure  1).  Details  of  the  burner  construc¬ 
tion  and  operation  are  given  elsewhere  (McEnally  and  Pfefferle,  1997).  The  oxidizer  was 
air  (710  cm3/min),  while  the  fuel  was  a  mixture  containing  methane  (240  cm  /mm)  an 
argon  (2.4  cm3 /min).  The  added  argon  aided  species  measurements  as  discussed  below. 
All  flowrates  were  governed  by  electronic  mass  flow  controllers  and  are  estimated  to  be 

accurate  to  within  5%. 

3.2  Temperature  and  Soot  Volume  Fraction  Measurements 

Temperatures  were  measured  with  uncoated  75  /im  wire-diameter  type  R  thermocou¬ 
ples.  By  measuring  profiles  in  coflowing  flames  with  various  size  thermocouples,  we  have 
determined  that  while  conduction  errors  of  up  to  100K  occur  at  some  flame  locations  with 
125/xm  wires,  they  are  negligible  with  75  /zm  wires  (McEnally  et  ai.,  1997).  The  measured 
junction  temperatures  were  converted  to  gas  temperatures  with  a  heat  balance  (radiation 
energy  loss  =  convection  energy  gain)  and  the  following  parameters:  Nusselt  number  - 
2  26,  junction  diameter  =  240  fim  (measured  with  a  microscope),  temperature-dependent 
thermal  conductivity  data  for  air  (Holman,  1986),  and  temperature-dependent  type  S 
thermocouple  emissivities  (Bradley,  1961).  We  estimate  that  the  absolute  error  m  the  gas 
temperatures  is  less  than  50K  in  soot-free  flame  regions.  In  soot-contammg  flame  regions, 
particles  rapidly  deposit  onto  the  junction,  rendering  the  junction  diameter  and  emissivity 
values  listed  above  incorrect.  Therefore  the  junction  was  cleaned  before  every  measure¬ 
ment  by  retracting  it  into  the  oxidizing  layer  outside  the  flame,  and  the  temperature  was 
recorded  quickly  (0.75  seconds)  after  the  junction  was  reinserted  into  the  flame. 


9 


Soot  volume  fractions  were  measured  with  a  new  technique,  thermocouple  particle 
densitometry  (TPD),  that  was  suggested  by  Eisner  and  Rosner  (1985)  and  recently  imple¬ 
mented  by  McEnally  et  al.  (1997).  A  clean  thermocouple  was  introduced  into  the  flame 
and  its  junction  temperature  was  recorded  periodically  for  200  seconds;  the  soot  volume 
fraction  was  then  determined  by  optimizing  the  agreement  between  the  measured  history 
and  one  calculated  with  the  thermophoretic  mass  transfer  formulation  of  Eisner  and  Ros¬ 
ner  (1985)  In  the  work  cited  above  we  demonstrated  that  TPD  measurements  agreed 
with  soot  volume  fractions  measured  by  the  conventional  technique  of  laser  extinction  in 
counterflow  and  coflowing  nonpremixed  flames.  We  developed  TPD  and  emp  oye  i  ere 
because  the  combination  of  low  soot  volume  fractions  ( <  1  ppm)  and  short  pathlengt  s  ( 
mm)  make  laser  extinction  measurements  difficult  in  methane  coflowing  flames  (Garo  et 
al.,  1986),  (Shaddix  et  al.,  1994).  With  TPD,  soot  volume  fractions  can  be  easily  measured 
down  to  at  least  0.02  ppm  by  allowing  sufficient  time  to  pass  for  a  measurable  change  in  the 
junction  temperature  to  occur.  The  technique  also  allows  spatially-resolved  measurements 
without  tomographic  reconstruction.  The  results  have  an  absolute  accuracy  of  ±  50%  and 
a  relative  accuracy  of  db  10%. 

Recent  experiments  have  shown  that  nonpremixed  flames  contain  two  kinds  of  soot 
particles:  an  early  type  that  appears  translucent  when  viewed  under  an  electron  micro¬ 
scope,  and  a  later  type  that  appears  dark  (Mergaridis  and  Dobbins  1989).  Similarly,  wo 
types  of  particles  have  been  identified  in  premixed  flames,  where  the  early  particles  were 
shown  to  be  transparent  to  visible  light  (d'Alessio  et  al.,  1992).  Results  from  our  labo¬ 
ratory  suggest  that  soot  volume  fractions  measured  by  TPD  include  both  particle  types, 
white  laser  extinction  and  laser-induced  incandescence  respond  to  only  the  dark  particles 
(McEnally  and  Pfefferle,  1996),  (Kdylii  et  al.,  1996). 

3.3  Species  Measurements 

Species  concentrations  were  measured  by  extracting  gases  from  the  flames  with  a 
narrow-tip  quartz  microprobe,  and  analyzing  them  with  on-line  mass  spectrometry.  The 
sampling  procedure  and  its  performance  are  discussed  m  greater  detai  m  McEnally  an 
Pfefferle  (1997)  where  it  was  concluded  that  the  microprobe  spatial  resolution  was  roughly 
one  millimeter,  reactions  in  the  microprobe  tip  have  only  a  minor  affect  on  the  measured 
concentrations,  and  condensation/adsorption  losses  to  the  sample  manifold  walls  were  in¬ 
significant  for  biphenyl  (C12H10;  154  AMU)  and  all  lighter  species.  Following  Pun  (1992), 
we  have  equipped  the  probe  with  an  oscillating  sapphire  fiber  that  extends  through  its 
orifice  and  grinds  down  any  soot  deposits  that  accumulate  there.  This  prevents  complete 
clogging  of  the  orifice,  enabling  us  to  make  measurements  in  particle-containing  flame 

regions. 

Most  species  were  quantified  with  a  custom-built  photoionization/time-of-flight  mass 
spectrometer  (PTMS)  that  can  simultaneously  measure  all  C3  and  larger  hydrocarbons 
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(except  propane  and  butane)  with  part-per-million  detection  limits  (Bermudez  and l  Pf- 
efferle,  1995).  This  instrument  generates  ions  via  single-photon  photoionization  wi  a 
low-intensity  vacuum-ultraviolet  (118.2  nm)  laser  beam.  Benzene  concentrations  were  di¬ 
rectly  calibrated  and  are  accurate  to  within  a  factor  of  two. 

A  commercial  electron-impact/quadrupole  mass  spectrometer  (EQMS;  Leybold  Inficon 
Transpector  C100M)  was  used  to  measure  methane,  acetylene,  and  argon.  The  EQMb 
ionizes  target  gases  with  a  fixed-energy  70  eV  electron  gun,  which  produces  extensive  ion 
fragmentation.  Methane  was  measured  with  its  CH+  fragment-ion  to  avoid  interferences 
from  the  0+  fragment-ions  of  water  and  carbon  dioxide,  argon  with  its  Ar  ion  to  avoid 
interferences  from  C3H+  ions,  and  acetylene  with  its  parent  ion.  The  methane  and  acety  ene 
concentrations  were  directly  calibrated  and  are  accurate  to  within  50%.  However,  ethylene 
fragment  ions  may  account  for  up  to  30%  of  the  acetylene  signal;  ethylene  and  its  C2H3 
fragment-ion  were  completely  obscured  by  nitrogen  at  all  flame  locations,  so  only  an  upper 
bound  of  2000  ppm  could  be  determined  for  its  mole  fraction. 

The  flowrate  into  the  microprobe  varied  with  flame  position  due  to  changes  m  the  gas 
temperature  and  partial  clogging  of  the  orifice.  To  account  for  this,  all  PTMS  and  EQMS 
mass  spectra  were  normalized  with  the  EQMS  argon  signal  by  assuming  that  the  argon 
concentration  was  uniform  throughout  the  flame.  As  mentioned  above,  argon  was  added  to 
the  fuel  so  that  its  concentration  was  1%  in  both  reactant  streams.  At  a  height  close  to  the 
burner  surface,  where  soot  clogging  was  insignificant,  the  radial  temperature  profile  could 
be  recovered  from  the  raw  argon  signals,  indicating  that  the  constant- argon  assumption  is 

valid. 

4.  Numerical  Method 

Our  goal  is  to  obtain  a  discrete  solution  of  the  two-dimensional,  elliptic,  governing 
conservation  equations  for  the  temperature,  gas-phase  and  soot  mass  fractions,  velocities 
and  vorticity  on  a  two-dimensional  mesh  the  initial  nodes  of  which  are  formed  by  the 
intersection  of  the  lines  of  the  mesh  Mr 

Mr  =  {0  =  r0  <  rx  <  . . .  <  ri  •  .  •  <  rMr  =  Ro},  C4-1) 

and  the  mesh  Mz 

Mt  =  {0  =  z0  <  zx  <  ...  <  zj ...  <  zM.  =  Z}.  (4-2) 

Computationally,  we  combine  a  steady-state  and  a  time-dependent  solution  method.  A 
time-dependent  approach  is  used  to  help  obtain  a  converged  numerical  solution  on  an 
initial  coarse  grid.  Grid  points  are  then  inserted  adaptively  and  Newton’s  method  is 
used  to  complete  the  problem.  Specifically,  the  spatial  operators  in  the  governing  partial 
differential  equations  are  approximated  by  finite  difference  expressions.  Diffusion  terms 
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are  approximated  by  centered  differences  and  convective  terms  by  upwind  approximations. 
With  the  difference  equations  written  in  residual  form,  we  seek  the  solution  U  of  the 
system  of  nonlinear  equations 

F(U)  =  0.  "  t4-3) 

For  an  initial  solution  estimate  U°  which  is  sufficiently  close  to  U' ,  the  system  of  nonlinear 
equations  in  (4.3)  can  be  solved  by  Newton’s  method.  This  leads  to  the  iteration 

J(Un)(Un+1  -  Un)  =  -\nF(Un),  n  —  0,1,2, -  (4-4) 

where  J(Un)  =  dF(Un)/dU  is  the  Jacobian  matrix  and  An  (0  <  A  <  1)  is  a  damping 
parameter. 

We  point  out  that  with  the  spatial  discretizations  used  in  forming  (4.3),  the  Jacobian 
matrix  in  (4.4)  can  be  written  in  block-nine-diagonal  form.  We  form  several  columns  of  the 
Jacobian  simultaneously  using  vector  function  evaluations  and  the  Jacobian’s  given  sparsity 
structure.  The  Newton  equations  are  solved  by  a  preconditioned  BiCGSTAB  itera  ion. 
The  Newton  iteration  continues  until  the  size  of  ||U"+1  -  U"h  »  reduced  appropriately 
Grid  points  of  the  two-dimensional  mesh  are  determined  by  equidistnbutmg  positive  weig 
functions  over  mesh  intervals  in  both  the  r  and  z  directions.  During  the  pseudo-transient 
portion  of  the  calculation,  the  size  of  the  time  steps  is  chosen  by  monitoring  the  local 
truncation  error  of  the  time  discretization  process.  Due  to  the  cost  of  forming  the  Jacobian 
matrices  with  detailed  transport  and  finite  rate  chemical  kinetics,  a  modified  Newton 
method  is  implemented  along  with  several  theoretical  estimates  that  determine  when  a  new 
Jacobian  should  be  reformed.  These  theoretical  results  help  increase  the  overall  efficiency 
of  the  algorithm  (see  Smooke  et  a i.,  1990,  1996  for  more  details). 

5.  Results  and  Discussion 

In  this  section  we  apply  the  model,  solution  procedure,  and  experimental  methods 
discussed  in  Sections  2,  3  and  4  to  a  confined,  axisymmetric,  coflow,  methane-air  diffusion 
flame  with  the  goal  of  helping  to  interpret  flame  structure  and  soot  formation  phenom¬ 
ena.  Initial  modeling  results  were  compared  with  the  experimental  values  for  bulk  flame 
parameters;  i.e.,  flame  height,  peak  flame  temperatures,  etc.  Initially,  flame  height  was 
overpredicted  by  the  model  by  15%  and  temperatures  were  underpredicted  by  about ,100- 
150K.  This  lack  of  agreement  was  considered  unacceptable.  Upon  analysis  of  both  the 
model  and  the  experiment,  it  became  apparent  that  conduction/radiation  from  the  base  o 
the  flame  to  the  burner  exit  tube  in  the  experiment  caused  noticeable  preheating  of  ot 
inlet  fuel  and  air  streams.  Calculations  demonstrated  that  the  flame  height  and  local  tem¬ 
peratures  were  significantly  affected  by  these  assumed  inlet  temperatures.  Increasing  the 
inlet  temperatures  from  300  to  420K  (mass  flow  held  constant)  provided  better  agreement 
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to  the  flame  height  and  local  temperatures  (see  later  discussion).  This  approach,  while 
not  perfect,  offers  a  tractable  computational  solution  to  a  challenging  problem.  This  in¬ 
crement  is  consistent  with  attempts  to  measure  the  temperature  at  the  burner  exit  plane. 
The  computational  results  demonstrate  the  importance  of  accurate  specification  of  ini¬ 
tial  conditions,  specifically  temperature  (and  resultant  effects  on  local  velocities),  on  the 
numerical  solutions. 

Time-stepping  from  nearly  converged  solutions  led  to  the  build  up  of  large  soot  par¬ 
ticles  at  the  base  of  the  flame  in  the  wake  of  the  burner  lip.  Neither  surface  growth  nor 
inception  can  produce  soot  in  this  low  temperature  region.  Several  numerical  experiments 
demonstrated  that  the  source  of  this  build-up  was  thermophoresis.  Small  soot  particles 
generated  a  little  further  from  the  lip  (at  higher  temperatures  and  inside  the  flame  front) 
were  driven  inside  this  small  wake  due  to  large  temperature  gradients.  Coalescence  of  these 
small  particulates,  continually  fed  into  the  wake  region  over  long  computation  times,  led  to 
the  formation  of  large  particles  which  were  not  flushed  away  by  the  flow.  In  heavier  sooting 
co-flow  flames,  such  build-up  with  time  has  been  observed  experimentally  as  evidenced  by 
the  occasional  requirement  to  shut  down  and  clean  particulates  from  the  exit  tube  (San¬ 
toro,  1996).  This  build-up  is  not  observed  experimentally  for  the  lightly  sooting  methane 
flame.  This  process  has  been  suppressed  in  the  model  by  setting  the  thermophoretic  ve¬ 
locity  to  zero  below  750K.  This  empirical  correction  does  not  cause  any  noticeable  change 
of  predicted  soot  levels  in  the  remainder  of  the  flow  field. 

The  chemical  kinetic  mechanism  used  in  this  work  was  derived  from  GRIMech  2.11  with 
its  recommended  thermodynamics  (Bowman  et  ai.,  1995).  Modifications  of  the  mechanism 
included  removal  and  addition  of  reaction  paths  which  are  believed  to  have  a  minimal 
impact  on  the  bulk  chemistry  in  the  system.  All  reactions  related  to  NO  formation  and 
other  odd  nitrogen  species  were  stripped  from  the  set.  In  addition,  a  series  of  reactions 
related  to  the  formation  and  oxidation  of  benzene,  and  related  species,  were  included 
(see  Table  1  of  Hall  et  ai.,  1997).  The  dominant  ring  forming  reaction  was  propargyl 
recombination.  This  is  consistent  with  material  balance  analyses  of  benzene  formation 
mechanisms  based  on  the  experimental  species  concentrations.  The  burner  fuel  tube  had 
an  inner  radius  of  Ri— 0.5556  cm  with  a  wall  thickness  of  0.0794  cm.  The  radius  of  the 
coflow  was  R0  =  4.7625  cm.  Conditions  at  the  fuel  and  oxidizer  jet  inlets  (plug  flow)  were 

such  that 

vF  =  5.52cm/sec,  vox  =  12.54cm/sec,  (5.1) 

Tf  =  420K,  Tox  =  420K.  (5.2) 

All  calculations  were  performed  on  an  IBM  RS/6000  Model  590  computer.  In  the  com¬ 
putations  to  be  presented  12  soot  size  classes  were  included  in  the  model.  Starting  from 
a  converged  solution  for  a  methane-air  flame  without  the  sectional  equations,  we  typi¬ 
cally  obtained  converged  solutions  for  the  complete  gas-soot  problem  in  several  hours  of 
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computer  time. 

Temperature  Predictions 

Examination  of  the  curves  in  Figure  2  demonstrates  that  the  model  predicts  the  general 
shape  and  structure  of  the  flame  (location  of  peak  flame  temperatures,  etc.)  reasonably 
well.  Temperature  is  perhaps  the  most  important  of  the  flame  parameters  as  temperature 
has  a  direct  influence  on  virtually  all  the  other  flame  properties  (e.g.,  flame  height  and 
width,  species  concentrations,  soot,  etc.)  A  more  careful  comparison  of  the  profiles  m 
Figure  2  shows  that  the  model  underpredicts  the  experimental  measurements  along  the 
centerline  by  about  100  K  at  all  heights  within  the  flame.  In  addition,  while  the  location 
and  value  of  the  peak  temperature  are  reasonably  well  described  up  to  about  a  height  of  z 
=  2  cm,  temperatures  are  overpredicted  at  the  outside  edge  (i.e.,  predicted  temperatures 
do  not  fall  off  as  rapidly  as  the  data  with  increasing  radius).  Finally,  the  expenmenta 
flame  closes  to  the  centerline  (at  z  «  40  mm)  more  rapidly  than  does  the  computed  flame. 
The  computed  low  centerline  temperatures  and  the  delayed  contraction  of  the  flame  tip, 
presumably,  are  very  much  related  to  the  inability  to  match  accurately  the  flame  height. 
These  deficiencies  in  the  model  are  attributable,  at  least  partly,  to  the  uncertainties  m  t  e 
inlet  boundary  conditions.  As  argued  below,  these  limitations  of  the  computed  temperature 
field  influence  the  prediction  of  gas  species  and  particulates.  One  possible  explanation  of 
the  low  temperatures  along  the  centerline  is  the  lack  of  inclusion  in  the  model  of  absorption 
by  methane  of  energy  radiated  from  the  flame  front. 

Prediction  of  Species 

Figure  3  shows  comparison  of  the  methane  profiles  as  a  function  of  radius  at  several 
heights  above  the  burner.  Again  the  general  features  of  the  radial  profiles  are  reproduced 
well  by  the  model.  Typically,  methane  is  overpredicted  by  a  factor  of  two  along  t  e 
centerline.  This  result  is  consistent  with  lower  temperatures  and  a  taller  computed  flame. 
Alternatively,  the  model  overpredicts  the  fall-off  with  increasing  radius,  most  noticeably  at 
the  lowest  height,  z  =  1.0  cm.  This  accentuated  fall-off  is  undoubtedly  a  direct  result  of  the 
higher  predicted  temperatures  in  the  wings  of  the  flame,  and  hence  more  rapid  chemistry. 
The  factor  of  three  overprediction  at  z  =  2.5  cm  on  the  centerline,  is  a  consequence  of  the 
(experimental)  rapid  decrease  in  methane  concentration  in  this  region  of  the  flame,  i.e., 
the  flame  tip,  and  the  inability  to  match  the  flame  height  exactly. 

Similar  comments  can  be  made  in  regard  to  the  acetylene  and  benzene  profiles  (Fig¬ 
ures  4  and  5).  In  these  cases,  however,  centerline  values  are  underpredicted,  presumably 
because  methane  decomposition  is  suppressed  at  the  low  predicted  temperatures  along  the 
centerline,  as  discussed  above.  Not  surprisingly,  the  relatively  small  underprediction  of 
acetylene  (factor  of  two)  along  the  centerline,  leads  to  a  larger  deficit,  m  the  computed 
benzene  profiles  (a  factor  of  three)  on  the  centerline.  Since  these  species  are  key  to  the 
growth  and  inception  of  particulates,  the  combination  of  their  low  values  as  well  as  the 
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low  temperature  predictions  along  the  centerline  must  lead  commensurately  to  low  pre¬ 
dictions  of  soot  in  these  regions.  Of  importance,  however,  is  that  the  concentrations  of 
these  species,  as  well  as  temperatures,  are  better  described  in  the  flame  front  region  most 
critical  to  inception  and  growth.  Until  a  satisfactory  solution  is  found  to-  predict  more  ac¬ 
curately  centerline  temperatures,  quantitative  predictions  of  formation  of  the  translucent 
droplets/particles  along  the  centerline,  such  as  observed  by  Megaridis  and  Dobbins  (1989) 
will  likely  not  be  achieved. 

Soot  Particulates 

Figure  6  depicts  the  computed  and  measured  values  of  the  soot  volume  fraction  (cm 
soot/cm3  gas)  as  a  function  of  radius  at  «  =  2.0,  2.25,  2.5  and  2.75  cm.  As  with  the 
species  profiles,  the  computations  describe  the  general  features  of  the  soot  field.  Peak 
soot  volume  fractions  are  underpredicted  by  a  factor  of  about  three,  but  the  shape  of 
the  soot  volume  fraction  profile  and  its  shift  towards  the  centerline  with  increasing  z  is 
reproduced  well.  The  experimental  peak  values  collapse  to  the  centerline  more  rapidly, 
as  discussed  above,  because  of  the  difference  in  the  computed  flame  height.  Centerline 
values  are  significantly  underpredicted  (order  of  magnitude)  because  of  the  combination  of 
predicted  low  temperatures  (100K),  low  acetylene  (factor  of  two)  and  low  benzene  (factor 
of  three).  It  is  reasonable  to  assume  that  should  these  bulk  flame  properties  be  more 
accurately  modeled,  our  ability  to  attain  quantitative  predictions  of  the  soot  field  would 
be  enhanced  dramatically. 

Predicted  soot  volume  fractions  are  strongly  dependent  on  predicted  acetylene  pro¬ 
files,  which  generally  are  dependent  upon  the  kinetics  supplied  by  GRIMech,  but  also  are 
depressed  when  scrubbing  is  included  in  the  calculations.  Without  the  added  chemical 
mechanisms  and  soot  formation  steps,  predicted  levels  of  acetylene  are  very  close  to  the 
experimental  values,  depending  on  the  flame  location.  Measured  peak  acetylene  concentra¬ 
tions  are  about  12,000  ppm.  The  measured  peak  soot  volume  fractions  are  4.7xl0-7  which 
is  approximately  equivalent  to  4000  ppm  acetylene.  Hence,  suppression  of  the  acetylene 
peak  by  the  soot  formation  is  expected  and  is  observed  in  these  calculations  -  even  in 
these  lightly  sooting  flames.  The  reduced  acetylene  levels,  in  turn,  suppress  the  levels  of 
soot  predictions,  since  both  the  inception  and  the  growth  processes  are  strongly  depen¬ 
dent  on  the  acetylene  concentration.  Predicted  acetylene  concentrations  more  consistent 
with  experimental  results  would  result  in  more  favorable  agreement  between  predicted  and 
measured  soot  volume  fraction  in  Figure  6.  Soot  volume  fraction  predictions  based  on  this 
growth  model  have  some  residual  sensitivity  to  the  number  of  size  classes  at  the  value  of 
12  employed,  which  was  the  largest  number  permissible  given  the  number  of  species  and 
grid  nodes  relative  to  available  computer  memory.  While  size  class  convergence  will  have 
to  be  addressed  ultimately,  factors  such  as  inlet /boundary  conditions  and  flame  chemistry 
may  have  a  much  greater  influence  on  the  accuracy  of  the  solution  and  its  ability  to  repro¬ 
duce  bulk  flame  parameters  (specifically,  temperature).  As  will  be  seen,  calculations  with 
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another  soot  growth  model  show  much  more  rapid  convergence  with  the  number  of  size 
classes. 

2D  Representations 

In  Figures  7a-d  and  8a  we  illustrate  contour  plots  of  the  same  variables  plotted  m 
Figures  2-6.  The  two-dimensional,  color  contour  of  computed  temperatures  m  Figure  7a 
clearly  shows  the  peak  temperature  regions  in  the  flame  front.  The  highest  temperatures 
within  the  flame  are  not  at  the  flame  peak,  but  in  the  outer  annulus  regions.  Peak  height 
can  be  inferred  from  these  profiles  by  locating  the  point  at  which  the  temperature  peaks  on 
the  centerline,  for  example.  The  contour  plot  of  methane  shows  a  dramatic  loss  of  methane 
very  low  in  the  flame.  This  extremely  rapid  decay  is  principally  due  to  the  rapid  diffusion 
of  molecular  nitrogen  to  the  centerline,  rather  than  the  reaction  of  methane.  (In  fact  it 
is  exactly  this  rapid  diffusional  process  that  has  led  to  confusion  in  the  literature  when 
attempts  have  been  made  to  separate  dilution  from  temperature  effects  and  their  relative 
importance  to  soot  formation.).  Predicted  peak  acetylene  and  benzene  profiles  occur  on 
the  centerline  about  a  centimeter  below  the  flame  tip  near  z  =  2.8  cm.  Benzene  seems 
to  decay  more  rapidly  above  this  location  than  does  acetylene.  The  rapid  decay  in  e 
acetylene  appears  more  closely  aligned  with  the  flame  tip.  Alternatively,  the  earlier  decay 
in  the  benzene  is  linked  to  the  soot  inception  process,  which  attains  its  peak  rates  m  this 

region  of  the  flame. 

Soot  Growth  and  Inception 

Figures  8b-d  contain  contour  plots  of  the  net  soot  production,  soot  inception  and  soot 
oxidation,  respectively.  Comparison  of  the  net  soot  production  (Figure  8b)  with  inception 
(Figure  8c)  shows  that  the  growth  region  is  surface-  rather  than  inception-dominated. 
Higher  in  the  flame  net  growth  gives  way  to  oxidation  as  the  particles  burn  out  (Figure 
8d).  The  surface  growth  contours  are  a  convolution  of  the  temperature  and  acetylene 
contours  (Figures  7a  and  7c),  strongly  reflecting  the  annular  shape  of  the  temperature 
profile,  while  inception  is  an  image  of  the  benzene  contour  of  Figure  7d. 

Soot  Oxidation 

The  rates  of  oxidation  of  the  particulates  by  OH  are  shown  m  Figure  8d.  Maximum 
oxidation  rates  by  02  (not  shown  here)  are  about  a  factor  of  5  to  10  below  those  by  OH. 
Oxidation  by  02  occurs  in  two  separate  regions  of  the  flow.  Near  the  base  of  the  flame, 
02  diffuses  underneath  the  flame  and  contributes  to  soot  oxidation  on  the  fuel-rich  side 
of  the  flame,  but  these  oxidation  rates  are  low  due  to  the  limited  oxygen  supply.  02  also 
contributes  to  soot  oxidation  in  the  post  flame  region,  but  again  at  low  rates  since  OH 
has  already  consumed  the  majority  of  the  particulates  as  they  pass  through  the  flame 
front.  The  hydroxyl  radical  is  clearly  the  dominant  oxidizing  species  in  this  flame.  Peak 
oxidation  rates  occur  in  the  conventional  wing  tips  of  the  flame  and  high  oxidation  rates 
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continue  to  the  centerline.  These  results  (i.e.,  the  strong  dominance  of  OH  as  the  principal 
oxidizing  species)  are  wholly  consistent  with  those  of  Garo,  et  al.  (1990)  who  analyzed  t  e 
dynamics  of  a  methane- air  coflow  diffusion  air  and  more  recently  of  Kennedy,  et  al.  (19  ) 

in  their  detailed  modeling  of  an  ethylene  coflow  diffusion  flame.  Together,  the  present 
and  prior  results  cast  doubts  on  the  validity  of  models  which  neglect  the  role  of  oxidation 
by  OH  in  coflow  flames  (for  example,  Kaplan,  et  al.  (1996)  and  Sivanthanu  and  Gore 
(1994)).  The  cause  of  the  high  importance  of  OH  as  the  principal  oxidizing  species  is  its 
high  superequilibrium  level  in  the  flame  front.  (See  discussion  below.) 

In  this  study,  the  soot  field  has  been  modeled  using  a  constant  0.13  collision  frequency 
for  reaction  after  Neoh,  et  al.  (1981).  In  the  Kennedy,  et  al.  (1996)  study,  this  reaction 
efficiency  was  varied  linearly  from  0.05  low  in  the  flame  to  0.2  in  the  tip  flame  region, 
to  assist  in  matching  soot  profiles.  In  addition,  Garo,  et  al.  interpreted  their  data  to  be 
consistent  with  a  varying  collision  frequency  with  height  in  the  flame.  While  it  is  apparent 
that  such  methods  could  assist  matching  the  experimental  and  modeling  profiles  of  soot, 
we  have  not  yet  adopted  this  approach. 

Superequilibrium  OH 

Smooke  et  al.  (1992)  have  computed  superequilibrium  levels  of  hydroxyl  in  methane 
air  coflow  flames.  Santoro  and  coworkers  (1994)  have  measured  OH  concentrations  and 
temperatures  in  ethylene  flames  and  found  these  radical  concentrations  to  be  affected  by 
the  level  of  soot  in  the  flame.  They  argued  that  high  oxidization  rates  of  soot  decreased  OH 
concentrations.  Such  effects  were  not  noticeable  in  this  study,  presumably  due  to  the  much 
lower  concentrations  of  soot.  The  OH  measured  by  Puri,  Santoro  and  Smyth  (1994)  were 
described  to  be  at  superequilibrium  levels.  To  examine  levels  of  superequilibrium  in  the 
present  study,  calculated  OH  (referred  to  here  as  “flame  OH”  have  been  compared  to  the 
equilibrium  levels  of  OH  based  on  1)  the  local  temperature  (constant  temperature  calcula¬ 
tion),  2)  the  full  equilibrium  (constant  enthalpy  calculation),  and  3)  a  partial  equilibrium 
assumption  based  on 

O2  ■+■  H%  2,0 H 

and  the  local  flame  temperature.  The  ratio  of  this  flame  OH  to  the  equilibrium  levels 
are  plotted  in  Figure  9  as  a  function  of  the  radius  at  *  =  4.01  cm.  The  maximum  soot 
oxidation  rate  occurs  at  this  height  at  a  radius  of  about  0.31  cm.  At  this  radius,  the 
OH  attains  concentrations  a  factor  about  10  above  the  equilibrium  value,  thus  validating 
further  the  existence  of  large  OH  superequilibrium  concentrations  as  an  important  con¬ 
tributor  to  the  destruction/burnout  of  soot.  Chemical  systems  which  cannot  support  such 
high  superequilibrium  levels  may  be  more  prone  to  soot  emissions.  The  partial  equilibrium 
assumption  (assuming  knowledge  of  the  H2,  02  concentrations  and  temperature)  provides 
a  much  better  estimate  of  hydroxyl,  at  least  through  the  flame  front  region  where  soot  is 

being  oxidized. 
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Particle  Size 


Figure  lOa-d  shows  the  evolution  of  soot  particle  size  distribution  in  space  as  a  result  o 
surface  growth  and  coalescence.  The  sectional  mass  boundaries  have  been  placed  linearly 
on  a  logarithmic  scale.  Peak  soot  mass  concentration  occurs  in  size  class  nine,  whose 
mass-weighted  mean  diameter  is  about  48  nm,  which  should  be  considered  an  estimate.  A 
larger  number  of  sections  than  that  employed  here  would  be  required  to  obtain  an  accurate 
average  particle  size.  This  average  size  is  clearly  larger  than  typical  determinations  o  - 
30  nm  in  laboratory  flames  and  suggests  some  refinements  in  our  soot  model  are  still 
required.  The  fact  that  the  suspended  mass  in  the  highest  class  remains  relatively  small 
shows  that  the  upper  boundary  of  the  particle  size  domain  has  been  set  high  enough  to 
avoid  unphysical  accumulation  of  particles  in  the  highest  size  class.  The  smallest  size 
classes  are  observed  to  disappear  rapidly  in  the  post  flame  zone,  while  a  portion  of  the 
largest  particles  survive  to  the  end  of  the  computational  domain. 

Particle  size  distribution  within  the  flame  can  also  be  visualized  from  Figures  11a  and 
11b  in  which  the  mass  fractions  of  each  of  the  12  soot  mass  classes  (me)  are  plotted  as  a 
function  of  radial  distance  for  heights  of  1  and  3  centimeters.  Low  in  the  flame,  the  sixth 
soot  class  contains  the  peak  mass  fraction  and  at  3  cm  (just  before  the  onset  of  rapid  soot 
oidation  by  OH  inside  the  flame  front),  the  maximum  occurs  in  the  ninth  class.  At  this 
higher  position,  substantial  soot  is  found  also  near  the  centerline  of  the  flame.  However, 
the  average  particle  for  this  “centerline”  soot  is  much  smaller  than  that  just  inside  the 
flame  front,  a  result  qualitatively  consistent  with  the  results  by  Megandis  and  Dobbins 
(1989).  Since  benzene  concentrations  (Figure  7d)  and  hence  inception  rates  (see  Figure 
8c)  are  peaking  in  this  region  of  the  flow,  the  particles  depicted  on  the  centerline  are  like  y 
formed  along  or  near  to  the  centerline  rather  than  transported  via  thermophoresis  from 
the  flame  front  to  this  centerline  region. 

Radiation  Effects 

While  it  might  be  expected  that  thermal  radiation  effects  would  be  small  in  a  lightly 
sooting  flame  like  that  investigated  here,  simulations  in  which  the  radiative  loss  term,  Eq. 
2  21  is  set  equal  to  zero  show  that  radiation  effects  are  significant.  Without  radiation, 
predicted  peak  temperatures  are  122K  higher,  a  change  large  enough  to  affect  predicted 
soot  levels  and  species  concentrations.  Predicted  peak  volume  fractions  increase  by  a 
factor  of  three,  with  peak  acetylene  and  benezene  levels  increasing  by  45%  and  2-%, 
respectively.  The  significance  of  thermal  radiation  in  this  lightly  sooting  flame  arises  from 
residence  times  that  are  relatively  long  compared  to  those  for  counterflow  diffusion  flames, 
in  which  radiation  effects  are  significant  only  for  extremely  low  strain  rates  (Garcia,  et 
al.,  1996).  With  the  relatively  low  values  of  soot  concentration  and  the  fact  that  the  soot 
does  not  occupy  the  peak  temperature  regions,  the  thermal  radiation  m  the  experiment  is 
estimated  to  have  comparable  contributions  from  the  gas  (water  vapor  and  carbon  dioxi  e) 
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and  the  soot.  In  the  simulations,  where  the  predicted  soot  levels  are  lower,  the  radiation  is 
dominated  by  the  molecular  species.  The  physical  dimensions  of  this  flame  are  sufficiently 
small  that  self-absorption  or  optical  thickness  corrections  should  be  minimal. 

Variations  of  Inlet  Gas  Temperatures 

Early  modeling  efforts  with  the  inlet  gas  temperatures  at  298K  resulted  in  a  flame 
height  about  15%  higher  than  the  experimental  value.  In  addition,  predicted  peak  flame 
temperatures  were  about  125K  below  the  experimental  value.  Since  soot  volume  fractions 
are  a  strong  function  of  temperatures,  these  low  predicted  temperatures  substantially  lim¬ 
ited  our  ability  to  predict  accurately  the  soot  concentrations.  Local  temperatures  m  the 
unburned  gases  just  above  the  burner  lip  were  measured  to  be  much  higher  than  298K, 
presumably  due  to  radiation  preheating  and  conduction.  Accurate  treatment  of  this  prob¬ 
lem  would  require  knowledge  not  only  of  temperatures,  but  also  of  local  velocities  at  each 
location  in  the  burner  inlet.  To  1)  analyze  this  problem  in  a  simplified  manner,  2)  provide 
information  relative  to  a  sensitivity  study,  and  3)  elevate  the  peak  flame  temperatures 
to  values  closer  to  the  experimental  levels,  calculations  were  performed  for  elevated  tem¬ 
peratures  in  the  entire  fuel  flow,  the  entire  air  flow,  or  both  streams.  For  each  case,  the 
mass  flows  per  unit  area  were  held  constant  and  hence  velocities  increased  with  higher 

temperatures. 

Table  1  compares  the  peak  soot  volume  fractions  and  the  peak  flame  temperatures 
obtained  for  several  solutions  at  inlet  temperatures  ranging  from  298  to  500K.  An  ex 
animation  of  these  numbers  demonstrates  that  an  increase  in  the  air  temperature  has  a 
much  more  dramatic  effect  on  both  the  peak  flame  temperature  and  volume  fraction  than 
comparable  changes  in  the  fuel  temperature.  Presumably,  this  effect  is  due  to.  the  higher 
mass  weighted  enthalpy  content  of  air  under  stoichiometric  conditions.  The  higher  flame 
temperatures  which  arise  from  the  elevated  inlet  air  temperatures  lead  to  a  proportionally 
greater  increase  in  the  soot  generation. 

Comparisons  to  Fairweather,  et  a 1.  model 

Finally,  the  model  was  rerun  with  the  inception  and  surface  growth  model  of  Fair- 
weather  et  a 1.  (1992).  Soot  volume  fraction,  net  soot  production,  inception,  and  oxidation 
are  illustrated  in  Figures  12a-d.  In  this  model,  which  has  been  applied  to  ethylene-  and 
acetylene-fuelled  coflow  flames  (Kennedy,  1996;  Sivathanu  and  Gore,  1994),  both  inception 
and  surface  growth  are  proportional  to  acetylene  concentration.  Its  inception  and  surface 
growth  rates  are  used  here  in  the  polydisperse,  sectional  representation  rather  than  the 
monodisperse  representation  of  the  original  model.  Additionally,  oxidation  by  OH,  not 
included  in  Fairweather,  et  al.  (1992),  has  been  accounted  for  in  this  implementation. 
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TABLE  1 


Effect  of  Fuel  and  Oxidizer  Inlet  Temperatures  on 
Peak  Soot  Volume  Fraction  and  Peak  Flame  Temperature 

Tj/Tox 

fv 

Tmax 

298/298 

3.30  x  lO'8 

1942  K 

298/420 

8.71  x  10~8 

1976  K 

298/500 

1.43  x  10"7 

2000K 

420/298 

4.67  x  10"8 

1948  K 

500/298 

5.87  x  lO"8 

1954  K 

420/420 

1.14  x  10"7 

1984  K 

500/500 

1.56  x  10"7 

2005  K 

The  predicted  soot  concentrations  (Figure  12a)  are  seen  to  be  slightly  larger  than  those 
predicted  with  the  base  model.  Given  the  previously  discussed  uncertainties  associated 
with  the  chemical  mechanism  and  burner  boundary  conditions,  however,  no  judgements 
about  the  merits  of  this  model  vis-a-vis  the  base  model  are  warranted.  Inception  is  more 
significant  in  the  Fairweather,  et  al.  (1992)  model  (compare  Figures  12c  and  8c),  and  this 
seems  to  give  calculations  using  it  much  less  sensitivity  to  the  number  of  soot  sections 
employed.  Predicted  soot  volume  fractions  using  five  and  twelve  sections  were  not  signifi¬ 
cantly  different.  The  spatial  distribution  of  inception  in  the  two  models  is  very  different, 
as  Figures  8c  and  12c  show.  The  Fairweather,  et  al.  (1992)  model  has  an  inception  profile 
that  is  strongly  annular,  reflecting,  like  the  surface  growth,  the  convolution  of  temperature 
and  acetylene  profiles.  Oxidation  in  this  latter  model  is  also  dominated  by  OH. 

Summary  of  Differences  Betweeen  Model  and  Experiment 

Differences  between  our  predictions  and  the  experimental  data  include:  overprediction 
of  flame  height,  underprediction  of  peak  temperatures,  too  slow  decay  of  methane  along 
the  centerline  and  a  slow  decay  of  temperature  in  the  wings  of  the  flame.  In  addition, 
the  model  significantly  underpredicts  the  peak  soot  volume  fractions  in  the  flame  and 
tends  to  overpredict  particle  size.  Much  of  the  differences  can  be  directly  linked  to  the 
inadequate  prediction  of  the  temperature  field  and  the  flame  height.  Contrasting  these 
results  are  the  very  good  agreement  of  the  temperature  field  and  flame  height  obtained  m 
the  modeling  of  a  methane  coflow  flame  in  a  study  of  NO  formation  Smooke,  et  al.  (1997). 
Principal  differences  between  this  previous  study  and  that  of  the  present  investigation 
includes  dilution  of  the  fuel,  higher  flow  rates,  and  a  lifted  flame.  We  believe  the  most 
likely  of  these  parameters  affecting  uncertainties  in  the  present  study  is  the  inability  of 
our  code  to  properly  model  the  attachment  (i.e.,  close  proximity)  of  the  flame  to  the 
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burner  lip  and  all  of  its  attendant  effects.  Other  possible  contributors  to  this  error  include 
uncertainties  in  the  kinetics.  For  example,  using  this  kinetic  mechanism  (hydrocarbon 
portion  of  GRIMech  2.11)  the  extinction  limit  in  a  counterflow  methane- air  diffusion  flame 
is  significantly  overpredicted  (380  cm"1  experiments  vs.  predictions  of  620  cm"1).  We  note 
also  that  there  are  still  fundamental  uncertainties  in  soot  growth  modeling  such  as  the  rates 
of  inception,  surface  growth,  and  coalescence  of  very  small  particles.  For  example,  a  50% 
increase  in  the  surface  growth  rate  constant  results  in  a  predicted  soot  volume  fraction  in 
much  better  agreement  with  data.  Until  the  temperature  field  is  more  accurately  described, 
it  is  difficult  to  offer  quantitative  arguements  about  the  accuracy  (or  lack  thereof)  of  our 
soot  model.  Nevertheless,  the  qualitative  results  and  conclusions  should  still  be  valid. 


6.  Conclusions 

The  axisymmetric  laminar  diffusion  flame  provides  a  natural  environment  m  which  one 
can  investigate  the  interaction  of  soot  formation  with  detailed  gas-phase  chemistry  in  a 
multidimensional  system.  In  this  paper  we  have  compared  experimental  measurements  of 
the  temperature,  soot  volume  fraction  and  selected  species  from  an  axisymmetric,  lami¬ 
nar,  methane-air  diffusion  flame  with  computational  results  obtained  by  generalizing  the 
sectional  model  presented  in  Hall  et  al.  (1997)  to  a  coflow  model.  This  work  represents  the 
first  time  that  a  detailed  chemistry/complex  transport  flame  model  has  been  coupled  to  a 
polydisperse  soot  growth  analysis  in  which  the  scrubbing  of  gaseous  species  and  both  gas 
and  particle  radiation  are  included.  The  detailed  chemistry  includes  species  up  through 
benzene  with  a  steady-state  model  for  the  formation  of  polyaromatic  hydrocarbons  to 
simulate  soot  inception.  Unfortunately,  descrepancies  between  the  experimental  data  and 
the  computed  results  for  several  bulk  flame  parameters  (such  as  temperatures  and  flame 
height)  limited  the  ability  to  quantitatively  evaluate  the  overall  validity  of  our  model(s). 
Soot  volume  fractions  calculated  using  either  our  growth/inception  formulation  or  that 
by  Fairweather  et  al.  (1992)  underpredicted  the  experimental  values  by  about  a  factor  of 
three.  These  low  values  are  directly  attributable  to  low  predicted  temperatures,  leading  to 
low  computed  values  of  critical  growth  species  such  as  acetylene  and  benzene.  Inaccuracies 
in  the  predicted  temperature  fields  may  be  partially  due  to  inadequate  modeling  of  flame 
attachment  at  the  burner  lip. 

Despite  these  descrepancies,  several  key  conclusions  can  be  extracted  from  this  study. 
These  results  strongly  support  earlier  studies  indicating  the  importance  of  soot  oxidation 
by  hydroxyl  radicals  and  is  largely  due  to  super-equilibrium  levels  of  OH  that  are  about  a 
factor  of  ten  above  equilibrium.  This  work  also  confirms  the  importance  of  fully  coupling 
radiation  into  the  flow  field  solutions.  Radiation  loss  reduces  the  peak  temperatures  by  over 
100K.  Ignoring  radiation  leads  to  increases  in  concentrations  of  key  soot  growth/inception 
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species  contributing  to  a  factor  of  three  increase  in  predicted  soot  volume  fractions.  Par¬ 
ticle  sizes  are  shown  to  vary  significantly  with  location  in  the  flame.  Particle  diameters 
grows  much  more  slowly  along  the  centerline  of  the  flame,  consistent  with  experimental 
measurements.  Alternatively,  while  nearly  all  particles  are  consumed  through  the  flame 
front  in  this  ’non-sooting’  coflow  flame,  a  few  small  particles  with  large  diameter  (and 
low  surface  area)  escape  to  be  oxidized  in  the  post-flame  region.  Finally,  a  comparison  of 
two  significantly  different  inception/growth  models  has  led  to  good  agreement  between  the 
predicted  soot  (volume  fraction  and  spatial  distribution). 
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Figure  Captions 


Figure  1  Schematic  of  the  burner  configuration. 

Figure  2  Computed  and  measured  temperature  profiles  at  various  heights  m  the  flame. 
Figure  3  Computed  and  measured  methane  profiles  at  various  heights  m  the  flame. 
Figure  4  Computed  and  measured  acetylene  profiles  at  various  heights  m  the  flame. 
Figure  5  Computed  and  measured  benzene  profiles  at  various  heights  m  the  flame. 
Figure  6  Computed  and  measured  soot  volume  fraction  profiles  at  various  heights  in  the 
flame. 

Figure  7a  Temperature  isotherms  of  the  computational  model  as  a  function  of  the  radial 
and  axial  coordinates. 

Figure  7b  Methane  isopleths  (mole  fractions)  of  the  computational  model  as  a  function 
of  the  radial  and  axial  coordinates. 

Figure  7c  Acetylene  isopleths  (mole  fractions)  of  the  computational  model  as  a  function 
of  the  radial  and  axial  coordinates. 

Figure  7d  Benzene  isopleths  (mole  fractions)  of  the  computational  model  as  a  function 
of  the  radial  and  axial  coordinates. 

Figure  8a  Soot  volume  fraction  isopleths  of  the  computational  model  as  a  function  of  the 
radial  and  axial  coordinates. 

Figure  8b  Net  soot  production  isopleths  of  the  computational  model  as  a  function  of  the 
radial  and  axial  coordinates. 

Figure  8c  Soot  inception  due  to  acetylene  as  a  function  of  the  radial  and  axial  coordinates 
(rate  of  equation  2.1). 

Figure  8d  Soot  oxidation  due  to  OH  as  a  function  of  the  radial  and  axial  coordinates 
(larger  negative  values  imply  higher  oxidation  rates). 

Figure  9  Profiles  of  the  ratio  of  the  OH  mole  fraction  to  the  equilibrium  levels  as  a  function 
of  the  radial  coordinate  at  z  =  4.01  cm. 

Figure  lOa-d  Soot  mass  fraction  isopleths  for  size  classes  7-10.  The  sectional  mass  bound¬ 
aries  (min=127.0,  max=1.5xl010)  have  been  placed  linearly  on  a  logarithmic  scale.  Peak 
soot  mass  concentration  occurs  in  size  class  nine,  whose  mass- weighted  mean  diameter  is 
about  48  nm  (see  Figure  11  for  the  mass  boundaries  of  each  class). 
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Figure  11a  Contributions  of  the  soot  mass  classes  to  the  total  soot  volume  fraction  at 
z=1.0  cm. 

Figure  lib  Contributions  of  the  soot  mass  classes  to  the  total  soot  volume  fraction  at 
z=3.0  cm. 

Figure  12a  Soot  volume  fraction  isopleths  of  the  Fairweather  et  a 1.  (1992)  model  as  a 
function  of  the  radial  and  axial  coordinates. 

Figure  12b  Net  soot  production  isopleths  of  the  Fairweather  et  a L  (1992)  model  as  a 
function  of  the  radial  and  axial  coordinates. 

Figure  12c  Soot  inception  isopleths  of  the  Fairweather  et  a 1.  (1992)  model  as  a  function 

of  the  radial  and  axial  coordinates. 

Figure  12d  Soot  oxidation  due  to  OH  as  a  function  of  the  radial  and  axial  coordinates  for 
the  Fairweather  et  a 1.  (1992)  model  (larger  negative  values  imply  higher  oxidation  rates). 
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COMPUTATIONAL  AND  EXPERIMENTAL  STUDY  OF  SOOT  FORMATION 
IN  A  COFLOW,  LAMINAR  ETHYLENE  DIFFUSION  FLAME 


Abstract 

A  sooting  ethylene  coflow  diffusion  flame  has  been  studied  both  experimentally  and 
computationally.  The  fuel  is  diluted  with  nitrogen  and  the  flame  is  slightly ^lifted .to  min¬ 
imize  the  effects  of  the  burner.  Both  probe  (thermocouple  and  gas-sampling  techniques) 
and  optical  diagnostic  methods  (Rayleigh  scattering  and  laser-induced  mcmidescence)  ar 
used  to  measure  the  temperature,  gas  species  and  soot  volume  fractions^  A  detaile 
growth  model  in  which  the  equations  for  particle  production  are  coupled_to  the  flow 
gaseous  species  conservation  equations  has  been  used  to  investigate  soot  formation  m 
flame.  The  two-dimensional  system  couples  detailed  transport  and  finite  rate  chemi^ry 
in  the  gas  phase  with  the  aerosol  equations  in  the  sectional  representation  _Theformu- 
lation  includes  detailed  treatment  of  the  transport,  inception,  surface  growth,  oxidation 
and  coalescence  of  soot  particulates.  Effects  of  thermal  radiation  and  particle  scrubbing  of 
gas  phase  growth  and  oxidation  species  are  also  included.  Predictions  and  measuremen  s 
of  temperature,  soot  volume  fractions  and  selected  species  are  compared  over  a  range  of 
heights  and  as  a  function  of  radius.  The  formation  of  benzene  is  primarily  controll 
the  recombination  of  propargyl  radicals  and  benzene  production  rates  are  found  tc >  h 
the  rate  of  inception  as  well  as  the  net  rate  of  soot  growth.  The  model  predicted  soot 
volume  fractions  well  along  the  wings  of  the  flame  but  underpredicted  ^t  vdume  frac¬ 
tions  by  a  factor  of  four  along  the  centerline.  Oxidation  of  part.cn  ate  is  dominated 1  by 
reactions  with  hydroxvl  radicals  which  attain  levels  approximately  ten  times  higher  than 
calculated  equilibrium  levels.  Gas  cooling  effects  due  to  radiative  loss  are  shown  to  have  a 
very  significant  effect  on  predicted  temperatures. 
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COMPUTATIONAL  AND  EXPERIMENTAL  STUDY  OF  SOOT  FORMATION 
IN  A  COFLOW,  LAMINAR  ETHYLENE  DIFFUSION  FLAME 


Introduction 

Combustion-generated  soot  particulates  from  land-based  sources  pose  a  significant 
health  risk  and  are  the  subject  of  stringent  new  EPA  regulations.  Now  soot  emissions 
from  aircraft  face  the  likelihood  of  tightened  regulation.  Besides  regulatory  issues  so 
contributes  to  thermal  radiation  loads  on  combustor  liners  and  turbine  blades_  Sootern 
sions  enhance  contrail  formation  and  tactical  visibility  of  military  aircraft  Further  i 
paction  of  soot  on  low  observable  surfaces  can  compromise  the  radar  signature  of  aircr  . 
Quantitative  understanding  of  the  soot  growth  and  oxidation  mechanisms  and  the  abili  } 
to  model  accurately  these  processes  may  be  critical  to  the  development  of  strategies  to 

control  emissions. 

Despite  the  complexities  of  modeling  soot  formation  in  flames  using  detailed  chemistrj', 
the  linkage  of  soot  production  to  radiation  and  bulk  flame  properties  is  so  strong  that 
the  coupled  treatment  of  this  problem  is  becoming  a  necessity  for  quantitative  mode  mg 
of  flame  structure.  Such  modeling  for  a  generic  multidimensional  configuration  is  still 
beyond  our  current  computational  ability.  The  laminar  diffusion  flame  however,  provides 
an  environment  to  investigate  the  interaction  of  soot  formation  with  detailed  gas-p  ase 
chemistry  in  a  multidimensional  system. 

Recent  investigators  [1,2]  have  modeled  jet  diffusion  flames,  using  simplified,  monodis- 
perse  soot  formation  models  with  skeletal  kinetic  mechanisms.  We  recently  modified  the 
sectional  soot  formation  model  developed  in  [3.4]  for  incorporat ion  into. 
nar  axisymmetric  diffusion  flame  (cylindrical  fuel  jet  surrounded  by  a  coflowinp  oxid .  ) 

[5].  This  model  emplovs  a  velocity-vorticity  formulation  [6]  in  which  the  go\  erning  c 
vation  equations  are  solved  with  detailed  transport  and  finite  rate  chemistry  submodels  to 
predict  the  temperature,  species  mass  fractions  and  velocity  fields  as  functions  of  the t*o 
independent  coordinates.  When  this  model  was  applied  to  a  sooting  methane-aii  flam  , 
comparisons  between  the  model  and  experiments  were  reasonable,  yet  thema  cimg 
bulk  flame  properties  was  insufficient  to  enable  quantitative  comparison  of  the  calculated 
and  measured  soot  profiles.  In  contrast  to  these  results  was  the  nearly  perfect  agreement 
obtained  between  the  model  and  optical  diagnostics  of  temperature,  fuel  ^  0  for  a 
non-sooting,  methane  coflow  diffusion  flame  [7].  This  latter  blue  flame  was  diluted  and 
lifted  far  from  the  burner  inlet.  A  speculated  problem  m  the  sooting  flame  which  was 
attached  to  the  burner  was  uncertain  inlet  boundary  conditions  due  to  preheating  of  e 

fuel  and  air. 

The  purpose  of  the  present  investigation  is  to  study  a  coflow  diffusion  flame,  but 
partially  lifted,  to  minimize  effects  of  uncertain  inlet  conditions  and  to  compare  results 
from  calculations,  intrusive  diagnostics  and  non-intrusive  diagnostics  to  determine  e 
strengths  and  weaknesses  of  each  approach.  Through  this  study,  it  is  expected  that  fla 
types,  amenable  to  modeling,  and  correspondingly  appropriate  diagnostic  methods  can 
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be  recommended  for  future  studies.  Furthermore,  this  study  should  provide  additional 
information  on  the  nature  of  soot  formation  and  growth  m  coflow  d,ffus.on  flames. 


Problem  Formulation 

Soot  Modeling 

Soot  kinetics  are  modeled  as  coalescing,  solid  carbon  spheroids  undergoing  surface 
growth  in  the  free  molecule  limit.  The  particle  mass  range  of  interest  is  divided  into  sec¬ 
tions  [8]  and  an  equation  is  written  for  each  section  including  coalescence,  surface  growth, 
and  oxidation.  Sectional  analysis  makes  it  possible  to  obtain  the  particle  size  distnbu  10 
without  a-priori  assumptions  about  the  form  of  the  distribution.  For  the  smallest  section, 
an  inception  source  term  is  included.  The  transport  conservation  equation  for  each  section 
includes  thermophoresis,  an  effective  bin  diffusion  rate,  and  source  terms  for  Sa®'P 
scrubbing.  The  gas  and  soot  equations  are  additionally  coupled  through  non-a 
radiative  loss  in  the  optically-thin  approximation.  The  inception  model  employe 
is  based  on  an  estimate  of  the  formation  rate  of  two-  and  three-rmged  aromatic  species 
(naphthalene  and  phenanthrene).  and  is  a  function  of  local  acetylene,  benzene,  phenyl  a 
molecular  hydrogen  concentrations  [5].  The  contributions  from  the  inception  processes 
are  incorporated  in  the  first  sectional  bin,  whose  lower  mass  boundary  is  set  equal  to 
mass  of  the  smallest  inception  species.  In  the  sectional  representation  [8],  the  sectio 
mass  boundaries  vary  linearly  on  a  logarithmic  scale.  The  number  of  sections  required  f 
convergence  must  be  examined  for  each  problem  and  depends  on  the  relative  magnitudes 
of  surface  growth  and  inception.  Oxidation  of  soot  by  02  and  OH  is  treated  as  describe 
in  [4].  The  surface  growth  rate  is  based  upon  that  of  Harris  and  Werner  [9]  with .a ^  nom¬ 
inal  activation  energy  of  31.8  kcal/mole  as  suggested  by  Hura  and  Glassman  [10].  Vie 
empirically  adjusted  the  Harris  and  Weiner  rate  by  a  factor  of  two  as  in  [5j. 

Governing  Equations  and  Numerical  Method 

The  asymmetric  computational  model  employs  the  gas-phase  diffusion  flame  equa¬ 
tions  in  the  velocity- vorticity  formulation  [6]  with  the  sectional  approach  presented  m  [  ]. 
Buoyancy  is  included  in  the  model.  The  result  is  a  strongly  coupled  set  of  elliptic  partial 
differential  equations.  We  solve  for  the  radial  and  axial  velocities,  the  vorticity-.  the  tern 
perature,  the  gas  phase  species  and  the  particle  sectional  mass  fractions  The  sjstem 
closed  with  the  ideal  gas  law  and  appropriate  boundary  conditions  are  applied  on  each  side 
of  the  computational  domain.  Local  properties  are  evaluated  using  vectorized  and  highl 
optimized  transport  and  chemistry  libraries  [11],  The  sectional  thermophoretic  vdocities 
in  the  free  molecule  regime  are  given  in  [3]  as  are  the  sectional  diffusion  velocities  which  a  e 
written  with  a  mass  weighted  mean  diffusion  coefficient  for  each  size  class.  In  the  optically- 
thin  radiation  model  used  in  our  calculations  the  significant  radiating  species,  m  addition 
to  particulates,  are  tf20,  CO  and  C02.  Given  the  length  scales  of  the  flame  investigated, 
it  is  highlv  unlikelv  that  self-absorption  is  important.  Although  the  soot  volume  fraction 
reaches  near  ppm  levels,  the  narrowness  of  the  soot  shell  (1-2  mm)  will  mitigate  any  self 

absorption  effects. 

The  governing  conservation  equations  are  solved  on  a  two-dimensional  mesh  by  combin¬ 
ing  a  steady-state  and  a  time- dependent  solution  method  [7].  A  time- dependent  approach 
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is  used  to  help  obtain  a  converged  numerical  solution  on  an  initial  coarse  grid.  Grid  points 
are  then  inserted  adaptively  and  Newton’s  method  is  used  to  complete  the  problem. 


Experimental  Methods 

Atmospheric  pressure,  overventilated,  axisymmetric,  coflowing,  nonpremixed  laminar 
flames  were  generated  with  a  burner  in  which  the  fuel  flows  from  an  uncooled  4.0  mm 
inner  diameter  vertical  brass  tube  (wall  thickness  0.038  mm)  and  the  oxidizer  flows  from 
the  annular  region  between  this  tube  and  a  50  mm  diameter  concentric  tube  (see  Figure  1). 
The  oxidizer  was  air  while  the  fuel  was  a  mixture  containing  ethylene  and  nitrogen.  Fuel 
flowrates  were  governed  by  electronic  mass  flow  controllers  accurate  to  within  5%.  The 
same  burner  apparatus  was  used  for  all  the  experiments.  The  temperature  of  the  brass 
tube  for  this  slightly  lifted  flame  was  less  than  330  K. 

Probe  Measurements 

The  probe  measurement  procedures  have  been  described  previously  [12,13].  Gas  tem¬ 
peratures  were  measured  with  75  pm  wire-diameter  Type  R  thermocouples  and  corrected 
for  radiation  heat  transfer  effects  using  standard  techniques  [12].  A  rapid  insertion  pro¬ 
cedure  was  used  to  minimize  errors  due  to  soot  deposition  onto  the  thermocouple.  In 
soot-free  regions,  the  absolute  uncertainty  of  these  measurements  is  estimated  to  be  ± 

K  and  the  relative  uncertainty  to  be  ±10  K. 

Soot  concentrations  were  measured  with  the  same  thermocouples  using  thermocouple 
particle  densitometry  (TPD),  a  technique  in  which  soot  volume  fraction  is  inferred  from 
measured  rates  of  soot  particle  mass  transfer  to  the  thermocouple  junction  [12],  The  results 
have  a  relative  uncertainty  of  30%,  and  an  absolute  uncertainty  of  50%. 

Species  concentrations  were  measured  by  extracting  gas  samples  from  the  flames  with  a 
narrow-tipped  quartz  microprobe  and  analyzing  these  samples  with  on-line  mass  spectrom¬ 
etry  [13].  Acetylene  and  ethylene  were  quantified  with  an  Extrel  C50  variable-ionization- 
energy  electron-impact/quadrupole  mass  spectrometer,  and  C3  to  C12  hydrocarbons  with 
a  custom-built  photoionization/time-of-flight  mass  spectrometer.  Measurements  were  di¬ 
rectly  calibrated  and  have  an  absolute  uncertainty  of  30,%. 

Profiles  were  generated  bv  moving  the  burner  with  translation  stages.  The  axial  and 
radial  coordinates,  designated  r  and  r,  have  a  relative  uncertainty  of  ±0.2  mm  and  an 
absolute  uncertainty  of  ±0.5  mm. 

Laser  diagnostic  measurements 

Using  planar  laser  imaging,  we  obtain  two-dimensional  fields  of  temperature,  fuel  con¬ 
centration,  and  soot  volume  fraction  in  the  C,HJN,  flame.  The  temperature  field  is 
determined  using  the  two  scalar  approach  of  Starner  et  al.  [14]  and  included  the  measure¬ 
ment  of  Rayleigh  scattering  and  the  use  of  the  computed  fuel  concentration. 

The  soot  volume  fraction  field  is  determined  by  laser-induced  incandescence  (LII).  At 
sufficient  laser  intensities,  the  LII  signal  has  been  shown  to  be  directly  proportional  to  soot 
volume  fraction  [15].  The  probe  measurements  of  the  soot  volume  fraction  are  used  for 
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calibration. 

The  second  harmonic  of  an  Nd:YAG  laser  (532  nm)  is  focused  into  a  18.0  mm tall  ver- 
tical  sheet  over  the  center  of  the  burner.  The  incandescence  and  scattered  light .sooB^ 
perpendicular  to  the  laser  axis.  The  light  passes  through  an  appropriate  interference  fi 
and  then  is  focused  onto  an  intensified  CCD  camera. 

For  Rayleigh  scattering,  images  at  two  downstream  locations  are  acquired.  FOT  th'  fast 
set  of  images  the  laser  sheet  is  3  mm  off  the  surface  of  the  burner  A  532  nm 
filter  (10  nm  FWHM)  is  used  to  collect  the  Rayleigh  scattering.  In  the  region  from  . 
22.0  mm  downstream,  interference  from  LII  and  particle  scattering  dominate  the  Rayleigh 
signal  which  is  not  plotted  in  this  region  (see  Figure  2).  Rayleigh  images  «  ^ 
with  the  bottom  edge  of  the  laser  sheet  just  above  the  flame  tip  (22  mm .  off  the  burne 
surface),  where  temperatures  are  just  below  adiabatic  flame  temperatures,  and  where  there 
are  no  interferences  from  LII.  Laser  energy  is  set  to  100  mJ/pulse. 

With  the  laser  sheet  S  mm  off  the  burner  surface,  laser-induced  incandescence  images 

are  acquired.  A  450  nm  (10  nm  FWHM)  interference  filter  is  used  to  collect  the  incan¬ 
descence  signal  In  the  region  of  greatest  incandescence  signal  a  survey  is  conducted  o 
incandes cence  signal  versuflaser  intensity  to  maximize  signal  intensity  without  saturation 
(e.g.,  soot  destruction)  at  any  point  in  the  flow  field.  A  1  images  are  corrected  for  optical 
throughput,  background  scattering  signals  and  non-uniformities  in  beam  profile.  Images 
are  also  corrected  for  flame  luminosity  and  non-uniform  detector  response. 


Results 

The  chemical  kinetic  mechanism  for  ethylene  combustion  has  45  species  and  233  reac¬ 
tions.  It  was  derived  from  GRIMech  1.2  [16].  based  upon  comparisons  to  experimental  data 
on  ethvlene  from  perfectly  stirred  and  flow  reactors  and  ignition  delay  data.  It  includes 
reactions  describing  the  formation  and  oxidation  of  benzene,  and  related  species. 

Fuel  and  nitrogen  are  introduced  through  the  center  tube  (4mm  id)  utilizing  a  parabolic 
velocity  profile  and  air  through  the  outer  coflow  with  a  plug  flow  profile.  Both  %e loci 
profiles  were  those  employed  in  the  experiments.  The  mass  fractions  at  the  burner  exi  axe 
0  32  and  0.6S  for  ethvlene  and  nitrogen,  respectively  and  the  bulk  averaged  velocity  is  3o 
cm/sec.  The  coflow  air  velocity  was  35  cm/sec.  Reactant  temperatures  were  assumed  to 
be  298  K.  All  radial  velocities  were  assigned  to  zero  at  the  flame  base. 

Calculations  were  performed  on  an  IBM  RS/6000  Model  590  computer.  In  the  com¬ 
putations  presented,  nine  soot  size  classes  were  included  m  the  model  with  approxim  ^ 
10,000  adaptively  refined  grid  points.  Starting  from  a  converged  solution  for  an  ethj 
air  flame  without  the  sectional  equations,  we  typically  obtained  converged  solutions  for  the 
complete  gas-soot  problem  in  several  hours  of  computer  time.  The  number  o  soot  bins 
in  these  calculations  was  constrained  by  the  maximum  memory  of  our  computer.  Based 
upon  the  relative  magnitudes  of  the  inception  and  surface  growth  rates,  we  anticipate 
the  restricted  number  of  bins  has  not  caused  significant  numerical  error. 

In  Figure  2.  temperatures  determined  from  the  model,  the  thermocouple  and  Rayleigh 
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scattering  are  compared.  Radial  comparisons  between  the  computations  and  thermocou¬ 
ple  measurements  at  several  axial  locations  are  plotted  in  Figure  3.  Agreement  between 
the  temperature  computations  and  measurements  is  excellent  throughout  the  flame.  The 
peak  temperatures  were  1953  K  for  the  computations  and  2040  K  for  the  thermocouple 
measurements.  Peak  Rayleigh  temperatures  were  somewhat  lower  than  the  computations. 
With  radiation  from  soot  suppressed,  the  peak  predicted  temperature  is  1990  K;  wltbbo 
gas  band  and  soot  radiation  suppressed,  the  peak  predicted  temperature  rises  to  2061  h. 
Integration  of  the  computed  radiative  dissipation  over  the  flame  volume  yields  a  predicte 
radiative  power  that  is  about  12%  of  the  total  heat  release.  Flame  height,  as  estimated 
by  the  attainment  of  the  peak  temperature  on  the  centerline  was  2.7  (±  0.05)  cm  tor  a 
three  cases.  This  agreement  is  dramatically  better  than  that  obtained  in  the  attached 
methane-air  diffusion  flame  [5]  in  which  the  experimental  inlet  conditions  were  not  well 
defined.  It  also  should  be  noted  that  the  excellent  agreement  on  the  air-side  of  the  flame 
contrasts  with  the  results  obtained  in  [5]  where  there  was  substantial  difference  between 
the  measured  and  computed  temperatures.  The  calculated  rise  m  temperature  along  e 
centerline  is  delayed  relative  to  the  increase  as  determined  from  the  probe  measurements, 
but  is  in  good  agreement  with  the  optical  results.  We  speculate  that  accurate  determina¬ 
tion  of  temperatures  in  this  region  of  the  flame,  where  thermal  gradients  are  very  large, 
with  a  thermocouple  is  difficult  due  to  conduction  along  the  thermocouple  wires. 


Ethylene  profiles  were  obtained  from  all  three  methods;  two-dimensional  contour  plots 
and  the  radial  profiles  at  various  axial  heights  (not  shown  here)  depict  excellent .  agree¬ 
ment  among  the  model  and  the  Raman  scattering  experiments.  Peak  concentrations  ot 
acetylene,  the  principal  carbon-containing  species  involved  in  surface  growth,  as  calcula  e 
by  the  model  and  as  determined  by  the  probe  measurements  are  4.3 7o  and  3.7%.  Peak 
concentrations  are  located  along  the  centerline  about  15  mm  above  the  burner  lip  an 
generally  the  profile  shapes  agree  well,  at  least  in  regions  of  the  flame  which  do  not  have 
other  contributions  to  the  mass  26  peak.  Experimentally,  acetylene  decays  a  little  faster 
along  the  centerline  than  does  the  model.  Contours  for  benzene  are  shown  m  Figure  4_ 
Peak  predicted  benzene  mole  fractions  are  l.SxHT4  versus  a  peak  experimental  value  of 
1.6xl0~4.  A  reaction  path  analysis  demonstrated  that  benzene  was  formed  principally 
through  propargyl  recombination  reactions  (to  form  H  +  C6H5).  Propargyl.  m  turn,  is 
formed  through  the  oxidation  of  diacetylene,  i.e.,  C4H2  +  OH  -*  C3H0  +  HCO  followe 
by.  C3H2  recombination  with  H  atoms.  When  soot  inception  scrubbing  of  be^ze^  1S _Jot 
activated  in  the  model,  the  peak  benzene  concentration  increases  by  a  factor  of  38  7o.  ibis 
contrast  demonstrates  the  importance  of  scrubbing  effects  and  demonstrates  how  benzene 
levels  are  determined  based  upon  a  balance  of  its  formation  and  its  conversion  to  higher 
molecular  weight  hydrocarbon  species.  Other  than  the  inception  process,  the  principal  loss 
mechanism  for  benzene  in  these  fuel-rich,  pyrolytic  zones,  is  via  the  thermal  decomposition 
of  phenyl  radical  and  such  decomposition  does  not  occur  rapidly  until  the  gas  temperature 
exceeds  1800  K,  a  point  after  the  bulk  of  soot  inception  and  growth  has  occurred. 


Soot  volume  fractions  /v,  as  determined  from  the  model,  thermophoretic  sampling, 
and  LII  (calibrated  based  on  probe  measurements)  are  illustrated  in  Figure  5.  Agreement 
between  the  two  experimental  techniques  is  considered  very  good.  Peak  soot  fv  from  t  e 
model  agrees  well  with  the  experimental  values;  however,  soot  still  peaks  off  the  centerline 
in  the  wing  region  for  the  model  and  on  the  centerline  for  the  two  sets  of  experiments.  Peak 


values  are  7.85  x  10'’  for  the  model  and  1.0  x  10-«  for  the  thermophoretic  measurements. 
When  the  predicted  values  ultimately  peak  along  the  centerline,  they  axe  a  factor  of  four 
below  the  measured  value.  As  illustrated,  the  thermophoretic  sampling  method  detects 
soot  particles  closer  to  the  burner  than  does  LII.  The  sampling  method  has  he  abflift 
to  detect  not  only  carbonized  soot  but  any  translucent  particles  as  well.  We  also  note  n 
Figure  5  that  the  computations  illustrate  an  extended  wing  region  compared  to  that  of  bo  h 
experimental  methods.  This  is  due,  in  part,  to  the  fact  that  the  soot  volume  fraction1 
the  lower  wing  regions  are  below  the  lower  detectivity  limits  of  the  experimental  methods. 
In  Figure  6  we  compare  computational  and  thermophoretic  sampling  soot  volume  factions 
as  a  function  of  the  radial  coordinate  at  several  axial  heights.  While  the  agreement  at  1 
mm  is  excellent,  it  is  clear  that  the  model  does  fall  short  in  being  able  to  predict  soot 
values  that  are  as  large  as  those  measured  on  the  centerhne.  The  tendency  f?r  so°* 
peak  in  the  wings  is  typically  observed  in  more  heavily  sooting  flames  [2].  Rckhvc  spatial 
distributions  for  the  separate  processes  of  surface  growth,  soot  inception  and  oxidation  as 
determined  from  the  model  are  illustrated  in  Figures  7,  8,  and  9,  respectively. 

A  brief  analysis  of  these  results  indicates  that  the  high  soot  oxidation  rates  observable  in 
the  wings  of  this  flame  are  attributable  to  super-equilibrium  OH  concentrat consistent 
with  our  previous  study  [5]  and  that  of  several  previous  investigations  [U]  The  lack  of 
including  the  effects  of  superequilibrium  OH  in  such  coflow  flames  are  likely  to  resu 
in  significant  errors  in  the  analysis/interpretations.  Variations  m  the  base  soot  model 
surface  growth,  inception,  and  oxidation  rates  were  carried  out  to  gain  understanding  of 
the  possible  causes  of  the  difference  between  the  model  and  the  experiments.  These  studies 
did  not  yield  a  clean  explanation  of  the  discrepancies.  We  conclude  that  the  ability  to  make 
quantitative  soot  predictions  remains  limited  by  some  fundamental  uncertainties  m  the  so 
model  (including  the  lack  of  aging  and  aggregate  formation  effects),  by  the  a  i  'J 
chemical  kinetic  mechanism  to  predict  accurately  the  concentrations  of  important  speoes 
(benzene,  propargyl,  acetylene  and  diacetylene)  and  possibly  by  the  lack  of  quantitat 
information  concerning  the  production  of  translucent  particles  [IS]. 


Conclusions 

A  slightly  lifted  ethvlene  jet  diffusion  flame  was  investigated  by  comparing  results 
from  two  sets  of  experimental  diagnostics,  one  of  which  was  intrusive  and  the  othernori- 
intrusive,  and  results  from  a  detailed  model  with  fully  coupled  equations  treating  rad 
and  soot  formation.  The  current  work  is  the  first  to  apply  a  detailed  chemistry  model  i 
a  multiple  section  soot  growth  model  to  a  flame  that  has  well  represented  burner/inlet 
conditions.  Agreement  among  the  experiments  and  computations  is  generalh  good  an 
some  cases  excellent.  In  addition,  by  comparing  the  results,  we  were  able  to  conclude  tha 
previously  identified  discrepancies  [5]  were  likely  the  result  of  uncertainties  in  the  burner 
inlet  conditions  when  the  flame  is.  attached  to  the  burner  lip.  Specifically  uncertain  inlet 
conditions  previously  led  to  over  prediction  of  the  flame  height  and  high  temperatures  m 
the  wings  of  the  flame.  For  the  lifted  flame,  the  model  was  able  to  reproduce  bulk  flame 
parameters  extremely  well,  including  flame  height,  species  concentrations  and  loca  tem¬ 
peratures,  given  some  uncertainties  in  the  experiments.  The  coupled  ^J^***^ 
in  this  and  in  previous  studies  reproduced  peak  soot  volume  fractions  to  within  -0%  but 
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had  some  difficulty  in  reproducing  accurately  the  distribution  of  soot  formed  along  the 
centerline  of  the  fkme  versus  that  formed  along  the  wings.  The  formation  of  benzene  as  a 
limit  to  the  inception  process  was  confirmed  in  this  study.  Benzene  formation  was  found  to 
be  governed  by  propargvl  recombination  and  propargvl  formation,  m  turn  was  controlled 

by  reactions  involving  diacetvlene. 
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Figure  Captions 

Figure  1.  Schematic  of  the  burner  configuration. 

Figure  2.  Temperature  isotherms  (K)  for  the  model  (left),  thermocouple  (center)  and 
Rayleigh  scattering  (right)  measurements. 

Figure  3.  Comparison  between  the  experimental  thermocouple  temperature  measure¬ 
ments  and  the  computed  temperatures  as  a  function  of  the  radial  coordma  e  a  sever 
axial  heights  (n=5  mm,  o=10  mm,  o=15  mm,  +=20  mm,  x-25mm). 

Figure  4.  Comparison  between  the  computed  and  experimental  benzene  mole  fraction 
isopleths. 

Figure  5.  Comparison  of  computed  (left)  and  experimental  soot  volume  fraction  isopleths. 
The  center  picture  contains  the  thermophoretic  sampling  measurements.  The  rig  gure 
contains  the  laser-induced  incandescence  measurements. 

Figure  6.  Comparison  between  the  computed  and  experimental  (thermophoretic  sampling 
technique)  soot  volume  fraction  as  a  function  of  the  radial  coordinate  at  several  axial 
heights  (o=18  mm,  +=20  mm,  x=22mm). 

Figure  7.  Computed  soot  surface  growth  isopleths  for  the  ethylene  diffusion  flame. 
Figure  8.  Computed  soot  inception  isopleths  for  the  ethylene  diffusion  flame. 

Figure  9.  Computed  soot  oxidation  (OH  +  0,)  isopleths  for  the  ethylene  diffusion  flame. 
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